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Final  Summary  Report 
on 

RESEARCH  AND  DEVELOPMENT  OF  HIGH 
TEMPERATURE  SEMICONDUCTOR  DEVICES 


FOREWORD 


The  principal  objectives  of  the  programs  conducted  under  Contracts 
NObsr  77532  and  NObsr  85424  were  to  carry  out  studies  and  to  investigate 
methods  and  techniques  which  would  lead  to  the  development  of  gallium 
arsenide  transistors.  This  work  was  to  include  the  following: 

(a)  Investigation  of  contact  materials,  etching  and  surface 
preparations,  Junction  forming  processes,  and  encapsulation 
techniques. 

(b)  Fabricate  feasibility  samples  of  a  high-temperature  low- 
power  gallium  arsenide  transistor. 

(c)  Fabricate  feasibility  samples  of  a  high  temperature  medium 
power  gallium  arsenide  transistor. 
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Final  Summary  Report 


on 

RESEARCH  AND  DEVELOPMENT  OF  HIGH 
TEMPERATURE  SEMICONDUCTOR  DEVICES 


ABSTRACT 


This  Final  Summary  Report  on  Department  of  the  Navy^  Bureau  of  Ships^ 
Contracts  NObsr  77532  and  NObsr  85424,  titled  Research  and  Development  of 
High  Temperature  Semiconductor  Devices,  covers  the  period  I  March  1959 
through  July  1962. 

Much  improvement  was  achieved  in  GaAs  single  crystals  used  for  device 

studies.  Single  crystal  GaAs  was  grown  by  the  Teal-Little  method  from 

ingots  compounded  by  the  gradient  freeze  technique.  Tunnel  diodes  were 

prepared  from  Zn  doped  degenerate  p-type  GaAs,  Transistors  were  generally 

fabricated  from  1  x  10*^  atoms/cm^  Sn  doped  n-type  material,  with  typical 

mobilities  in  the  range  of  3000  to  4000  cm^/vol t-sec.  Improved  crystal 

growing  techniques  yielded  material  with  dislocation  densities  as  low  as 
2 

3000-4000/cm  ,  an  order  of  magnitude  better  than  generally  obtained.  Suc¬ 
cessful  attempts  were  made  to  make  a  grown  diffused  GaAs  transistor  with 
current  gains,  3  ^  0.2. 

Initially  alloyed  contact  studies  were  made  to  obtain  both  ohmic  and 
rectifying  contacts  on  GaAs.  Ultimately,  studies  were  made  to  find  a  high 
efficiency,  high  temperature  emitter  alloy  for  the  transistor.  Au-Sn 
and  Au-Se  were  satisfactory  as  ohmic  contacts  on  p-  and  n-type  GaAs, 
respectively.  Only  Sn  and  Au-Sn,  of  several  alloys  investigated  as 
rectifying  contacts,  exhibited  satisfactory  injection  efficiency;  neither 
had  the  necessary  high  temperature  capability  (up  to  400®C)  for  the  transistor. 
Both  ohmic  and  rectifying  contacts  were  incorporated  and  alloyed  successfully 
on  the  transistor,  with  the  upper  temperature  limitation  of  the  unit  being 
restricted  to  300**C  or  less. 
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Zn>  Cd,  Mn,  S,  Mg;  and  Hg  diffusion  in  GaAs  for  forming  p-n  junctions 
was  studied.  Zn  diffusion  was  investigated  in  some  detail  because  Zn 
behaves  in  an  anomalous  manner.  The  diffusion  coefficient  of  Zn  varies  by 
about  4  to  5  orders  of  magnitude  over  a  relatively  small  concentration 
range.  Cd  diffusion  was  also  quite  interesting  because  of  its  apparent 
diffusion  in  a  direction  opposite  to  its  concentration  gradient.  We  have 
shown  this  is  primarily  because  of  the  independent  diffusion  of  the  decay 
daughter  of  the  radioactive  Cd  used.  Cd  like  Zn  seemed  to  be  sharply 
concentration  dependent.  Mn  diffusion  data  utilizing  electrical  methods 
are  reported.  The  diffusion  coefficients  of  Mn  are  considerably  higher 
than  those  of  Zn  at  equivalent  concentrations  and  appear  to  obey  a  co¬ 
error  function  distribution.  Sulphur  diffusion  data  using  both  electrical 
and  radioactive  tracer  results  are  shown.  Surface  compound  formation  in 
this  experiment  was  greatly  reduced.  Diffusion  coefficients  were  con¬ 
siderably  higher  than  reported  in  other  published  data.  Mg  results  as 
obtained  from  electrical  measurements  are  presented.  No  concentration 
dependence  of  the  diffusion  coefficients  was  observed;  but  Mg  obeyed  an 
error  function  distribution.  Hg  diffusion  data  were  obtained  by  tracer 
techniques  and  diffusion  was  found  extremely  slow  and  obeyed  a  complementary 
error  function. 

Zn  and  Mg  diffusion  in  the  presence  of  sputtered  SiO^  film  on  GaAs 
was  carried  out.  The  was  found  to  be  a  diffusion  mask;  depending  on 

film  thickness;  diffusion  time  and  diffusion  temperature.  It  also  pre¬ 
vents  decomposition  of  the  GaAs  surface  during  the  diffusion  cycle.  This 
film  may  lead  to  ultimate  planarization  of  GaAs  transistors. 

Thermal  conversion  of  GaAs  during  diffusion  cycles  was  definitely 
attributed  to  Cu  diffusion  and  means  of  controlling  the  degree  of  thermal 
conversion  was  found. 

Because  mechanical  polishing  damaged  GaAs  surfaces  and  caused  adverse 
effects  to  device  characteristics;  a  chemotaxial  polishing  technique  was 
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developed.  An  etch  consisting  of  80:10:10  parts  of  H^SO^,  and  H2O 

respectively  was  found  to  optically  polish  the  As  surfaces. 

Several  stains  and  etches  that  delineated  both  diffused  and  alloyed 
Junctions  in  GaAs  were  discovered;  e.  5HC1  and  IHNO^  in  a  3  to  2 
ratio  wi th  H^O. 

A  token  effort  of  post-alloy  diffusion  on  GaAs  with  a  Au-Sn-Zn  alloy 
yielded  transistors  with  a  3's  1. 

A  surface  *'passi vation*'  etch  of  SNaOH  to  IH2O2  was  observed  to  improve 
current  gains,  3,  by  as  much  as  an  order  of  magnitude.  Devices  niaintained 
this  value  of  3  when  exposed  to  ambients. 

Tunnel  diodes  were  fabricated  from  Zn  doped  p-type  degenerate  GaAs. 

Peak  current-to-val ley  current  ratios  as  high  as  70/I  were  realized  with 

voltage  swings  of  0,4  volts.  An  electrolytic  etch  was  evolved  that 

allowed  tailoring  the  diodes  to  desirable  peak  currents.  Peak  current 

2 

densities  of  5OO-5OOO  amperes/cm  were  observed,  with  capacitances  of 
typically  2-4  picofarads.  High  temperature  characterization  of  the  diodes, 
up  to  350^C,  was  successfully  carried  out,  and  and  values  vs  tempera¬ 
ture  presented.  Anomalous  humps  or  inflections  observed  in  the  I-V  curves 
of  some  units  were  attributed  to  deep-lying  impurity  levels  in  the  energy 
gap  into  which  a  second  tunneling  phenomena  could  have  been  taking  place. 

GaAs  transistors  were  proved  feasible.  N-P-N  mesa  units  that  exhibited 

current  gains  as  high  as  3OOO,  with  typical  3's  of  5^*20,  were  fabricated. 

- 1 0  -8 

Other  typical  dc  parameters  were  =10  -  10  amps;  »  15  volts; 

BVceo  “  ^5  volts;  and  =*  0.5  to  1.0  volt. 

Measurement  of  ’V  parameters  at  1000  cps  yielded  calculated  matched 
power  gain  values  as  high  as  43.6  db  in  the  common  emitter  case.  Typical 
hjb  50  ohms  with  h^j^  =  0.5  micromhos  were  measured. 
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parameter  measurements  vs  frequency  were  made  for  h.  and  . 

X  e  f  e 

The  h.  vs  frequency  curves  yielded  values  of  r'  C  measurements,  h. 
le  0  c  fe 

vs  frequency  yielded  values  for  the  gain  bandwidth  product  fj,  with  the 
highest  observed  being  800  Me.  Minority  carrier  lifetime,  in  the 

base,  was  calculated  from  the  f.  and  8  values  and  found  to  be  in  the  range 

-9  -8  ‘ 

of  10  ^  -  10  °  sec. 

Saturated  switching  time  measurements  were  performed  on  a  number  of 
units.  The  best  transistor  switched  10  ma  of  current  in  a  total  time  of 
10  nsec,  with  less  than  a  nanosecond  being  realized  in  the  storage  time. 

Beta  vs  temperature,  up  to  350®C,  was  observed  on  several  units.  On 
a  particular  unit  the  3  changed  from  14  to  11  over  the  temperature  range 
of  25°C  to  350®C.  No  permanent  degradation  in  the  unit  resulted  upon 
returning  to  room  temperature. 

Twenty-one  state-of-the-art  transistors  were  placed  on  storage  test 
at  150°C  for  1000  hours.  Again  no  serious  degradation  problem  was  observed 
in  any  of  the  units. 
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Final  Summary  Report 


On 

RESEARCH  AND  DEVELOFHENT  OF  HIGH 
TEMPERATURE  SEMICONDUCTOR  DEVICES 


I.  GALLIUM  ARSENIDE  MATERIALS  RESEARCH 


During  the  period  of  these  contracts,  considerable  improvement 
was  achieved  in  the  gallium  arsenide  single  crystals  used  for  device 
studies.  Although  this  materials  work  was  not  directly  supported  by 
contract  funds,  some  of  the  results  of  the  research  are  reported 
here  to  give  a  better  understanding  of  the  progr*ess  gained  in  transistor 
performance. 

Primarily,  the  gallium  arsenide  material  was  routinely  compounded 
by  the  gradient  freeze  technique,  which  in  many  cases  yielded  large 
single  crystal  areas.  Degenerate  n-  and  p-typa  ingots  were  prepared 
for  tunnel  diode  studies.  Some  single  crystal  areas  of  the  non¬ 
degenerate  n-type  ingots  were  used  in  technology  and  transistor 
fabrication  studies.  In  these  cases,  the  material  was  off  the  <11 1> 
orientation  and  needed  orientation.  In  general,  these  regions  were  of 

relatively  high  purity,  ^lO^^-lo'^  atoms/cc,  with  mobilities  in  excess 
2 

of  5000  cm  /volt  sec  at  room  temperatures. 

Mostly,  the  technology  and  device  fabrication  studies  were  per¬ 
formed  on  crystals  grown  by  the  Teal-Little  method.  They  were  pulled 
from  material  first  compounded  in  the  gradient  freeze  furnace.  Gen¬ 
erally,  these  crystals  were  grown  in  the  <11 1>  growth  direction.  See 
Fig.  1  for  a  picture  of  such  a  crystal.  These  crystals  were  grown 
both  non-degenerate  and  degenerate  and  in  intentionally  doped  and 
undoped  states. 


-I 
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Fig.  1  Gallium  Arsenide  Crystal,  III  Growth  Direction 


Table  I  shows  the  results  of  pulled  crystals  grown  from  melts 
doped  heavily  with  zinc  when  we  attempted  to  obtain  degenerate  p- 
type  GaAs  for  tunnel  diodes.  Despite  heavy  doping,  no  difficulty 
was  encountered  in  obtaining  single  crystals.  The  values  shown  in 
Table  I  were  obtained  from  Hall  measurements. 

Crystals  were  grown  to  study  the  effective  distribution  coeffi¬ 
cient  of  tin  in  gallium  arsenide;  (tin  is  an  n-type  dope  and  used  almost 
exclusively  in  crystals  for  transistor  fabrication). 

The  Hall  coefficients  near  the  tops  of  the  crystals  were  measured 
and  used  to  compute  the  donor  concentrations  in  the  solid.  Table  II 
shows  the  measured  values  and  the  resulting  segregation  coefficients. 

»3 

An  average  value  of  k  -  3.2  x  10  was  obtained.  The  crystal  result¬ 
ing  from  the  most  heavily  doped  melt  showed  tin  inclusions  which  in¬ 
dicated  that  the  solubility  limit  was  reached  at  a  value  of  about 
1  o 

7  X  10*  atoms/cc.  Experiments  on  tunnel  diodes  indicate  that  consider¬ 
ably  higher  tin  concentrations  are  obtained  by  alloying  at  lower  tempera¬ 
tures.  This  is  probably  because  retrograde  solubility  is  usually  more  pro¬ 
nounced  for  impurities  having  a  low  segregation  coefficient.  Mobility 

values  most  generally  realized  on  10^^  atoms/cc  tin  doped  crystals 

2 

ranged  from  3000  to  4000  cm  /volt-sec. 

Experiments  similar  to  those  described  for  tin  were  performed  to 
determine  the  segregation  coefficient  of  manganese  (found  to  be  a  p- 
type  impurity)  at  the  melting  point  of  gallium  arsenide.  Figure  2 
shows  the  dependence  of  the  segregation  coefficient,  k,  on  the  total 
manganese  concentration  added  to  the  liquid.  Figure  3  shms  k  as  a 
function  of  solid  concentration.  The  curves  show  that  k  varies  in¬ 
versely  as  the  square  root  of  the  solid  concentration  over  a  wide  range 
of  concentrations. 


-3- 


Table  I 


Degenerate  p-Type  Zn  Doped  GaAs  Crystals 


Crystal 

Carrier  Concentration 

Mobility 

Number 

(Cm"^) 

2 

(cm  /volt-sec) 

1 

2.6  X  10’^ 

79 

2 

8.2  X  10’^ 

58 

3 

1.2  X  10^° 

53 

Resistivity 

(ohm-cm) 

0.0058 

0.0025 

0.0018 
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Table  II 


Data  on  Tin->doped  Gallium  Arsenide  Crystals 

Melt  Composition  Excess  Donor  k  « 

(%  tin) _  Concentration  (cm  ) 

1.98 
7.85 
0.250 
0.113 
0.145 
0.107 


2. 12  X  1o’® 

6.9  X  lo’® 
1.97  X  1o’^ 
8.87  X  1o’® 
7.27  X  1o’® 

8.9  X  1o’® 


3.9  X  lO"^ 
3.0 

3.5 

3.0 

1.9 


10  10  10 

SEGREGATION  COEFFICIENT  k-^ 

Dependence  of  the  segregation  coefficient  of  manganese 
in  gallium  arsenide  on  the  Concentration  in  the  melt 
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Fig,  3  Dependence  of  the  segregation  coefficient  of  manganese  in  gallium 
arsenide  on  solid  concentration 


Crystal  growing  techniques  were  also  studied  for  their  effect  on 
dislocation  density  and  homogeneity  of  impurity  distribution.  Vari¬ 
ation  of  the  growth  parameters  -  pull-rate,  spin-rate,  doping  and 
crystal  diameter  -  resulted  in  improved  crystals  for  transistor  fab¬ 
rication.  Some  crystals  showed  dislocation  densities  as  low  as  3000  - 
4000/cm  ,  an  order  of  magnitude  lower  than  generally  obtained.  Im¬ 
purity  homogeneity  was  so  significantly  improved  that  flat  diffused 
Junctions  were  obtained  in  the  transistors. 


Attempts  were  also  made  to  apply  the  grown  diffused  technique  to 
gallium  arsenide  for  transistor  fabrication  with  tin  and  manganese  as 
Impurities.  The  general  procedure  was  first  to  grow  approximately 
half  the  crystal  from  a  charge  doped  with  a  small  amount  of  tin,  yield¬ 
ing  n-type  gallium  arsenide  with  properties  appropriate  for  the  tran¬ 
sistor  collector.  The  remainder  of  the  charge  was  re-etched  and  the 


impurities  added.  Amounts  calculated  to  give  solid  concentrations  of 

17  I S 

5  X  10 and  6  x  10  atoms/cc  of  manganese  and  tin,  respectively,  were 
used.  The  doped  charge  was  remelted^  and  the  initial  crystal  used  as 
the  seed  from  which  growth  was  continued  for  an  additional  2-10  minutes. 


Since  diffusion  studies  show  that  manganese  diffuses  much  faster  than 
tin,  manganese  should  form  a  p-type  base  region  ahead  of  the  heavily 
doped  n-type  emitter.  In  a  few  cases,  an  n-p-n  structure  was  actually 
formed  and  fabricated  into  transistors  with  current  gains,  ot,  as  high 
as  0.2.  Primarily^  this  technique  appeared  to  offer  difficulty  in 
obtaining  a  thin  base  region.  It  was  discontinued  because  of  later 
successes  with  the  n-p-n  mesa  structure. 
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II.  ALLOYED  CONTACT  STUDIES 


The  alloyed  contact  studies  were  initially  concerned  with  ob¬ 
taining  both  ohmic  and  rectifying  contacts  on  gallium  arsenide.  The 
ultimate  objective  was  to  find  a  high  efficiency^  high  temperature 
emitter  alloy  for  the  transistor. 

A,  Ohmic  Contacts 


Initially^  the  studies  were  to  find  low  resistance  ohmic  con¬ 
tacts  on  n-  and  p-type  gallium  arsenide.  Equipment  was  designed  and 
constructed  to  accurately  measure  the  specific  resistance  of  these 
contacts.  Basic  electrical  circuitry  and  apparatus  are  shown  in  Figs. 

4  through  8  and  are  described  in  the  following  paragraphs.  This  equip¬ 
ment  allowed  both  square-wave  ac  and  dc  iDeasurements . 

Direct  Current  Supply  (Fig.  4)--  A  set  of  cascaded  rheostats 
permits  continuous  adjustment  of  the  current  between  9  ma  and  3  e  (when 
a  1  ohm  sample  is  in  the  sample  holder).  A  switch  allows  the  source 
polarity  to  be  reversed. 

Square-wave  Alternating  Current  Supply  (Fig.  5)--  A  150  cps 
square-wave  was  available  from  a  Heathkit  S-3  electronic  switch.  This 
wave  has  less  than  0.  Ip^^ec  rise  time  with  negligible  ti1t»  but  the  out¬ 
put  is  a  high  impedance  voltage  source.  Therefore,  we  designed  and 
constructed  a  three-stage  transistor  current  amplifier  to  act  as  a 
current  source.  The  first  stage  serves  as  a  limiter  and  clamper.  The 
Input  square-wave  drives  the  transistor  between  1^^^^  and  the  saturation 
current.  The  second  stage  matches  the  impedance  between  the  limiter  and 
the  output  stage  and  prevents  the  latter  from  reaching  saturation.  Charge 
storage  is  prevented  and  distortion  reduced;  base  current  is  limited  by 
a  470  ohm  resistor^  which  is  bypassed  by  a  lOO^f  capacitor  to  prevent 
excessive  signal  transient  times.  The  output  stage  is  operated  between 


-9- 


6VD.C. 

STORAGE 

BATTERY 


Fig.  5  Power  amplifier  for  contact  resistance  studies 


an  output  current  determined  by  a  5  ohm  potential  divider  and  the  load 
resistance  and  I  .  Low  input  impedance  allows  a  minimum  of  collector 

CO 

current,  I  .  The  stage,  then,  is  simply  switched  on  and  off  by  the 

CO 

emitter  signal  of  the  second  stage.  A  SOOOiJif  capacitor  feeds  a  2  ohm 
load  and  allows  less  than  10  percent  tilt  of  the  wave,  but  prevents 
appreciable  dc  polarization  of  the  sample.  Up  to  I  ampere  of  square- 
wave  ac  to  a  2  ohm  load  is  supplied  by  current  amplifier.  Rise  times 
are  about  50M*sec  and  the  symmetry  is  excellent.  Distortion  causes 
less  than  0.2  percent  error  in  potential  measurements. 

Probe  and  Carriage  (Fig.  6)--  The  sample  is  placed  with  the 
end  contacts  bridged  to  nickel  blocks  by  indium  held  under  pressure 
between  the  jaws  of  a  lateral  vise.  Electric  current  contacts  are  made 
to  the  nickel  blocks.  The  vise  is  attached  to  a  movable  stage  which 
may  be  translated  laterally  in  0.01  mm  intervals  and  transversely  to 
align  the  sample  under  the  probe.  The  probe  is  an  0,040-in.  diameter 
tungsten  wire,  initially  sharpened  by  grinding,  and  finally  etched 
electrolytical ly  to  a  fine  point;  it  can  be  mechanically  lowered  to 
the  sample  and  held  in  place  by  a  small  adjustable  spring. 

Measurement  Selector  --  Two  poles  of  switch  A,  Fig.  7,  allow 
either  ac  or  dc  to  flow  through  a  1  ohm  precision  resistor  and  the 
sample.  Switch  B  permits  potential  measurements  to  be  made;  (1) 
across  the  precision  resistor  for  current  determination,  (2)  across 
the  sample,  (3)  between  the  front  sample  contact  and  probe,  or  (4) 
between  the  probe  and  rear  sample  contact.  Switch  B  is  wired  to  permit 
these  operations  without  reversing  the  K-3  input.  Switch  C  in  position 
“2**  allows  a  60-cycle  test  signal  to  pass  through  a  47  K  ohm  protective 
resistance  across  the  probe  contact.  When  probe  ohmic  contact  is  es¬ 
tablished,  the  oscilloscope  vertical  input  is  effectively  grounded  and 
a  horizontal  trace  displayed.  In  position  "1,"  the  switch  allows  the 
wave  form  of  the  potentials  to  be  monitored.  The  remaining  2  poles  of 
switch  A  permit  output  of  potentials  to  either  the  detector  circuit  or 
K-3  potentiometer. 
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Fig.  6  Photograph  of  the  probe  and  carriage 
for  contact  studies 


-13- 


Detector  Circuit  (Fig.  8)--  With  the  toggle  switch  in  peak- 
to-peak  position,  the  detector  is  a  conventional  voltage  doubler  cir¬ 
cuit,  loaded  by  a  1  M-ohm  load.  A  Pi  resistance-capacitance  network 
serves  as  a  ripple  filter  with  the  toggle  switch  in  peak  position;  the 
detector  is  a  shunt  rectifier  loaded  by  a  1  M-ohm  load.  Another  1 
M-ohm  resistor  isolates  the  Pi  filter  to  prevent  the  first  capacitor 
from  acting  as  a  voltage  divided  with  the  Pi  filter  capacitance.  In 
this  condition  the  detector  reads  potentials  resulting  from  current 
flowing  unidi rectional ly  through  the  sample.  Potentials  correspond¬ 
ing  to  reverse  current  flow  may  be  read  by  throwing  the  wafer  switch, 
which  reverses  the  diodes.  Thus,  though  square-wave  ac  flows  through 
the  sample,  direct  current  measurements  may  be  made. 

The  method  involves  passing  a  current  through  a  sample  of 
known  geometry  and  dimensions  and  determining  the  potential  at  various 
points  on  the  sample  by  precision  probing  techniques.  From  the  po¬ 
tential  drop  at  the  contacts,  the  sample  current,  and  the  cross- 
sectional  area,  a  value  of  contact  resistance  was  calculated.  Table 
III  shows  the  results  of  fourteen  different  contacting  materials  in¬ 
vestigated  initially  as  ohmic  contacts  on  n-  and  p-type  gallium  arsenide^ 
respectively.  It  should  be  noted  that  several  of  the  contacts  were  non- 
ohmic,  as  shown  by  differences  in  the  values  obtained  with  the  current 
flow  in  the  direction  from  the  contact  into  the  sample  from  those  ob¬ 
tained  when  the  current  flow  was  from  the  sample  to  the  contact.  Figure  9 
shows  a  typi ca I  poten ti a  I  profile  curve  f  rom  whi ch  the  contact  resi s  tance 
values  were  calculated.  Conclusions  from  these  and  subsequent  experiments 
showed  Au-Zn  and  Au-Se  to  be  satisfactory  ohmic  contacts  on  p-  and  n-type 
gallium  arsenide,  respectively.  They  are  now  used  in  the  fabrication 
of  state-of-the-art  transistors. 

Contact  resistances  determined  from  square-wave  ac  measure¬ 
ments  were  not  significantly  different  from  the  values  determined  from 
dc  measurements.  It  was  concluded  that  theomoelectric  effects  were 
negligible  in  GaAs.  The  ac  measurement  system  later  proved  useful  in 
studying  electrical  contacts  to  Bi^Te^  and  PbTe. 
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Detector  circuit  for  contact  resistance  studies 


Table  III 


Results  of  Contact  Specific  Resistance  Measurements 


2 

Specific  Resistance  (ohm-cm  ) 
GaAs  Current  into  Current  out 
Type  Contact  Sample  of  Sample 


n 

Ni  plate 

25.7 

0.440 

n 

Ni  plate  (preanodized) 

78.8 

17.1 

n 

Ni  plate  (sintered) 

0. 168 

0.184 

n 

Au  (alloyed) 

0.516 

0.516 

n 

Au-Sb 

95-5  (alloyed 
and  Ni  plated) 

0.250 

0.141 

n 

Au-Te 

95-5  (alloyed 
and  Ni  plated) 

0.058 

0.058 

n 

Au-Sb 

99.4-0.6  (al loyed» 
Ni  plated  and 
sintered) 

2. 

X 

lo"^ 

2. 

X 

10-3 

n 

Ag-Sn 

50-50  (alloyed) 

3.8 

X 

10-3 

3.8 

X 

10-3 

n 

Ag-Sn 

50-50  (alloyed 
and  Ni  plated) 

7.6 

X 

lo-'* 

7.7 

X 

lo-"* 

n 

Ni  plate  (sintered  and 
replated) 

6.36 

X 

10-3 

4.08 

X 

10-3 

n 

Ag-In 

99-1  (alloyed) 

6.28 

X 

lo”^ 

2.91 

X 

10-3 

n 

Au-In 

99-1  (alloyed) 

4.9 

X 

10-3 

3.9 

X 

10-3 

n 

Rhodium  plate  (sintered) 

0.50 

3.1 

X 

lO"^ 

P 

Ni  plate  (sintered) 

5.4 

X 

lo"** 

5.0 

X 

10“** 

P 

Ni  plate  (sintered 
and  replated) 

2.1 

X 

10-'* 

1.4 

X 

lo-'* 

P 

Au-2n 

99-1  (alloyed) 

5.0 

X 

lo"** 

1.6 

X 

10-3 

P 

Au-Zn 

99-1  (alloyed 
and  Ni  plated) 

6.8 

X 

10"3 

1.2 

X 

lO"^ 

P 

Rhodium  plate  (sintered) 

7.3 

X 

10-3 

2.1 

X 

lo"^ 

P 

Ag-AI 

72-28  (alloyed) 

7.5 

X 

10-3 

2.8 

X 

10-3 

P 

Ag-Mn 

95-5  (alloyed) 

2.7 

X 

10-^ 

1.6 

X 

10-^ 

P 

Ag-Mn 

95-5  (alloyed 
and  Ni  plated) 

1.9 

X 

10-^ 

9.5 

X 

10-^ 
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VOLTAGE 


Fig.  9  Potential  vs  distance  from  a  contact  made  by  alloying  i 
an  evaporated  film  of  95%  Au  -  5%  Sb  into  N-type  gallium 
arsenide  and  nickel  plating  surfaces 
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B.  Rectifying  Contacts 


Studies  of  rectifying  alloyed  contacts  on  gallium  arsenide 
were  designed  primarily  to  find  a  high  injection  efficiency^  high 
temperature  emitter  alloy  on  p-'type  material.  Many  alloys  such  as 
Au-Te»  Au-Sn«  Au-Se»  Ag-Sni  Sn-Te^  Sn«Se,  Si-Sn,  Au-Ge-Sb»  and  Sn 
were  investigated.  The  method  for  determining  the  high  injection 
efficiency  of  the  contacts  was-  first  to  observe  the  forward  I-V  trace 
of  the  diodes  on  a  curve  tracer^  then  attempt  to  fabricate  a  tran¬ 
sistor  incorporating  the  alloy  as  an  emitter. 

Pure  Sn  was  found  to  be  a  donor  and  a  good  standard  for  com¬ 
parison  with  other  alloys  for  injection  efficiency.  Diodes  were  fab¬ 
ricated  on  10^^  atoms/cc  p-type  gallium  arsenide  with  pure  Sn  and  were 
found  to  exhibit  the  characteristic  0.8  Volt  forward  breakover  voltage 
for  gallium  arsenide.  Transistors  made  with  Sn  emitters  exhibited 
current  gains,  8,  as  high  as  3000,  certainly  indicating  good  injection 
efficiency.  However^  Sn,  with  its  low  melting  point  of  232^C,  limited 
the  maximum  temperature  operation  of  the  transistors  far  below  the  ex¬ 
pected  350®C. 

Of  the  previously  mentioned  alloys,  the  only  other  one  that 
exhibited  desired  diode  and  transistor  characteristics  was  Au-Sn  62%  - 
38%  by  weight.  Diodes  made  with  this  material  had  the  same  character¬ 
istics  as  Sn  and  transistors  with  8's  of  about  100  were  fabricated. 
Again,  however,  the  maximum  operating  temperature  of  the  transistor 
was  limited  to  about  300^C.  The  62%Au  -  38%Sn  alloy  melts  at  4l8'’C, 
but  alloying  with  gallium  arsenide  lowers  this  melting  point  con¬ 
siderably. 


Search  for  an  ultimate  good  alloy  emitter  is  being  continued 
and  all  state-of-the-art  gallium  arsenide  transistors  are  now  being 
fabricated  using  the  Sn  emitters. 
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C.  Evaporated  Contacts 


The  goal  in  this  study  was  to  evaporate  and  alloy  all  the 
necessary  contacts  on  the  transistor.  The  ohmic  base  and  collector 
contacts  were  successful ly  evaporated  and  alloyed.  The  base  contact 
consisted  of  Au-Zn  and  the  collector  contact,  of  Au-Se,  Much  effort 
was  expended  on  evaporating  and  alloying  the  62%Au  -  38%Sn  alloy  as 
the  emitter  before  we  found  that  it  did  not  exhibit  the  desired  high 
temperature  properties. 

Initially,  we  experienced  difficulty  in  obtaining  good 
melting  and  alloying  of  the  Au-Sn  alloy  over  the  entire  area  of 
contact.  This  resulted  in  a  surface  barrier  contact  with  a  much 
lower  contact  potential  than  the  actual  p-n  junction,  with  which  it 
is  in  parallel,  and  therefore  yielded  0.5  volt  forward  breakover 
voltages  instead  of  the  characteristic  0. 8-1.0  volt  for  GaAs.  How¬ 
ever,  this  problem  was  finally  resolved  by  proper  preparation  of  the 
GaAs  surface,  just  prior  to  the  evaporation.  Stripes,  1  x  3  n^ils  and 
3  microns  thick,  were  evaporated  and  alloyed.  They  yielded  0.6  micron 
thick  regrowth  regions  over  the  entire  area  of  contact.  The  alloying 
cycle  was  carried  out  at  GOO^'C  in  a  reducing  atmosphere.  The  diodes 
obtained  exhibited  the  expected  0.8  to  1  volt  breakover  voltage  on  the 
forward  I-V  curves. 

It  was  not  until  this  study  was  completed  that  we  found  that 
the  62%Au  -  38%Sn  alloy  would  not  withstand  temperatures  in  excess  of 
300'’C  after  alloying  with  the  GaAs.  Work  will  be  continued  in  this 
area  when  a  suitable  high  temperature,  high  injection  efficiency  emitter 
alloy  is  found. 
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III.  DIFFUSION  STUDIES  IN  GALLIUM  ARSENIDE 


A.  General 


Diffusion  of  zinc^  cadium^  manganese^  sulfur^  magnesium,  and 

mercury  in  gallium  arsenide  have  been  studied  and  are  reported  here. 

Diffusion  of  zinc  has  been  studied  in  some  detail  because  zinc  seems 

to  behave  in  an  anomalous  manner.  The  diffusion  coefficient  of  zinc 

varies  by  some  4  to  5  orders  of  magnitude  over  a  relatively  small 

I  2 

range  in  concentration.  ' 

Cadmium  diffusion  is  also  quite  interesting,  because  it 
seems  to  diffuse  in  a  direction  opposite  to  its  concentration  gradient; 
that  is,  **uphi I l-di ffusion"  seems  to  occur.  We  have  shown  that  this 
results  primarily  from  the  independent  diffusion  of  the  decay  daughter 
of  the  radioactive  cadmium  used.  A  mechanism  that  seems  to  explain 
these  results  assumes  an  outdiffusion  of  dislocation  generated  vacan¬ 
cies  from  the  interior  of  the  sample  while  the  cadmium  is  diffusing 
inward. 


Manganese  diffusion  results  utilizing  electrical  methods 
are  also  reported.  The  diffusion  coefficients  of  manganese  are  con¬ 
siderably  higher  than  those  of  zinc  at  equivalent  concentrations. 

We  believe  this  is  because  manganese  diffuses  intersti tial ly,  although 
the  electrically  active  manganese  is  probably  substitutional;  this 
indicates  an  interstitial-vacancy  reaction. 

Sulfur  diffusion  results  using  both  electrical  and  radioactive 
tracer  results  are  also  shown.  Our  values  for  the  diffusion  coefficient 
are  much  higher  than  those  obtained  in  the  only  other  published  data.^ 

We  found  very  serious  surface  erosion  occurred  unless  suitable  pre¬ 
cautions  were  taken.  Our  data  were  taken  where  the  surface  erosion  was 
reduced  considerably,  and  we  believe  the  reported  values  are  accurate. 

Magnesium  results  are  presented  as  obtained  from  electrical 
measurements.  No  concentration  dependence  of  the  diffusion  coefficient 
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of  magnesium  was  observed,  but  it  obeyed  a  complementary  error  function 
distribution  in  GaAs.  No  tracer  data  could  be  obtained  on  magnesium 
to  substantiate  these  observations,  because  of  the  extremely  short  half- 
life  of  the  isotope. 


Mercury  diffusion  data  was  secured  by  the  tracer  technique 
and  found  to  obey  the  complementary  error  function  distribution  with 
extremely  low  values  of  diffusion  coefficients. 

In  the  theory  section  that  follows,  we  develop  a  broad  based 
diffusion  theory  from  first  principles.  With  this  theory,  based  on 
the  diffusion  of  an  atom  by  several  possible  modes  simultaneously,  we 
can  examine  the  following  special  cases:  (I)  an  atom  that  can  exist 
and  diffuse  in  more  than  one  state  of  ionization,  (2)  an  atom  diffusing 
by  both  interstitial  and  substitional  modes,  (3)  an  atom  diffusing 
primarily  by  an  interstitial  mode,  and  (4)  an  atom  diffusing  primarily 
by  a  substitutional  mode. 


B.  Theory  of  Diffusion  by  Parallel  Modes 

If  we  propose  that  diffusion  may  occur  by  several  modes  in 
a  crystal,  equations  can  be  developed  to  treat  the  situation.  We 
choose  to  treat  diffusion  at  a  given  region  of  the  crystal  as  being 
characterized  by  an  effective  diffusion  coefficient,  D,  and  we  assume  as 
a  first  approximation  that  Fick^s  first  law  is  obeyed.^  Hence  the  flux 
J,  of  a  given  type  of  atom  in  a  crystal  is  given  by: 


J 


(I) 


where  N  is  the  total  concentration  of  the  diffusing  atom,  and  x  the 
distance.  If  diffusion  occurs  by  several  modes,  the  flux  of  each 
mode  can  be  taken  into  account  with  the  expression 

dN. 

J  =  E  Jj  »  -E  Dj  ^  ,  (2) 

We  ignore  any  effect  from  interaction  of  concentration  gradients. 
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where  Jj>  Dj,  and  Nj  are  the  flux,  diffusion  coefficient,  and  concentra 
tion,  respectively,  of  atoms  in  the  i'th  mode. 


For  diffusion  by  two  parallel  modes,  this  becomes: 
dNj  aN2 


or,  remembering  that  N  =  N,  +  Nj,  we  have 

aN,  aN2 

afr  *^2  air^ 


M 

ax 


Now  the  term  In  brackets  can  be  identified  with  the  effective 
diffusion  coefficient  of  Eq.  (I),  so: 


D  =  D 


=  (0. 


aN, 

I  "W 


an 

aiJ 


^2 

1  .  n 


(3) 

{3a) 


If  the  fraction  in  each  mode  is  not  a  function  of  concentra- 
tion>  this  expression  becomes: 


N  N 

“  =  “l  N  ^ 


N 


=  (D,  ^  D^)  . 


(4) 

(4a) 


Equations  (3)  and  (4)  may  now  be  used  to  develop  several 
special  cases.  Equation  (4)  will  generally  be  used  because  of  its 
simplicity;  although  the  more  exact  expression;  (3);  is  also  calculated 
in  the  first  case  discussed. 


1 .  Diffusion  by  Different  Charge  States 

As  an  example  of  the  use  of  the  above  equations;  consider 
an  atom  in  a  semiconductor  that  can  be  either  un-ionized  (neutral);  or 
ionized  (charged).  The  fraction  of  the  atom  in  each  state  will;  in 
general;  be  a  function  of  the  impurity  concentration;  or  in  the  language 
of  band  theory;  the  fraction  is  a  function  of  the  position  of  the  Fermi 
level.  Further;  it  should  be  clear  that  in  general  the  diffusion  co¬ 
efficient  of  the  atoms  in  different  states  of  charge  will  be  different. 


-23- 


Following  is  a  discussion  of  the  diffusion  of  an  ac¬ 
ceptor  which  becomes  negatively  charged  when  ionized.  The  treatment 

of  a  donor  would  proceed  similarly.  The  acceptor  is  chosen  specif i- 
4 

cally  because  Allen  proposed  the  above  mechanism  to  explain  the  large 
concentration  dependence  of  the  diffusion  coefficient  of  zinc  in  GaAs. 
Unfortunately;  he  did  not  develop  his  equation  so  that  an  effective 
diffusion  coefficient  could  be  calculated.  We  show  in  the  following 
that  this  mechanism  predicts  a  diffusion  coefficient  variation  that 
is  in  very  poor  agreement  with  the  experimental  results  of  Cunnell 
and  Gooch^  and  with  our  work. 

The  reaction  by  which  an  un-ionized  acceptor,  A®, 
becomes  ionized  by  giving  up  a  hole,  e^,  (or  stated  equivalently, 
"accepting  an  electron  from  the  valence  band")  is 

A  A  +  e  . 

Reiss^  has  shown  that  if  the  acceptor  concentration  is  not  too  high 
this  type  reaction  can  be  treated  using  the  laws  of  chemical  equilibria, 
namely  the  Law  of  Mass  Action.  The  Mass  Action  equation  governing 
the  above  reaction  is 

N"  P 

_ _  L/  /.V 


where  end  ere  the  neutrel  end  ionized  ecceptor  concentretions^ 

P  the  hole  concentretion,  end  en  equilibrium  constent  given  by: 

Ny 

Ka  =  —  ®^P  ^-(Ea  -  Ey)/kT]  .  (6) 

Hy  Is  the  density  of  stetes  in  the  velence  bend,  Ea  the  ionization 
energy  of  the  ecceptor,  the  energy  et  the  top  of  the  valence  band, 
k  Boltzmann's  constant,  T  the  absolute  temperature,  and  the  factor  2 
takes  account  of  spin  degeneracy.  This  latter  factor  can  be  omitted 
if  desired  and  included  in  a  modified  ionization  energy. 
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Utilizing  two  other  relations,  namely  the  expression  for 
electron-hole  recombination  and  the  charge  neutrality  condition,  we  can 
calculate  the  fraction  in  each  state  at  any  acceptor  concentration* 

With  these  results  we  can  calculate  an  effective  diffusion  coefficient 
from  Eq*  (3)  or  (4).  The  relations  needed  to  complete  the  set  are  from 
the  electron-hole  equilibria. 


2 

np  =  ,  (7) 

where  n^  is  the  intrinsic  electron  concentration,  the  charge  neutrality 
condition  (ignoring  any  compensating  donors), 

n  +  N~  =  p  ,  (8) 

and  conservation  of  species 

^  •  (9) 


These  equations  can  be  manipulated  to  give  the  following  expression 
involving  N^o  and  N^: 


This  is  a  cubic  equation  in  N^e,  and  not  easily  solved,  but  at  very  low 
acceptor  concentrations,  i.e.,  «  n^  the  fraction  un-ionized,  as  well 
as  the  partial  derivative,  is 


e 

A 


n,  +  K, 


(n) 


Further,  at  moderate  and  high  acceptor  concentrations  where  »  n^ 
we  have  for  the  un-ionized  fraction: 
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2 


1/2 


2N.  +  K. 
A  A 

2N* 


and  for  the  partial  derivative 


.  (12) 


(13) 


The  effective  diffusion  coefficient  can  be  obtained  by 
inserting  (12)  or  (13)  into  (4)  or  (3)»  respectively.  We  have  done  this 
for  the  acceptor  zinc  in  GaAs  at  1000*^0  and  compared  it  with  the  experi¬ 
mental  results  of  Gunnell  and  Gooch.  ^  The  parameters  used  in  the  calcu¬ 
lation  were  n.  «  3.8  x  10^^  cm*"^,  ®  2.43  x  10^®  cm"^, 

*20-3  A  V  M 

and  N..  «  1.3  x  10  cm  •  These  values,  except  for  the  zinc  ionization 

*  ^ 
energy,  were  calculated  using  the  data  of  Ehrenreich.  The  choice  of 

zinc  ionization  energy  is  somewhat  arbitrary,  since  there  is  no  reliable 

"if 

published  value. 


The  results  of  the  calculation  are  shown  in  Fig.  10^  com¬ 
pared  with  the  experimental  results.  The  diffusion  coefficient  of  the 

4 

un-ionized  species  is  assumed  to  be  dominant.  Allen  claimed  qualita¬ 
tive  fit  with  this  theory  and  the  data  shown,  but  it  is  clear  that 
quantitatively  the  agreement  is  very  poor. 

2.  Parallel  Interstitial  and  Substitutional  Diffusion 

Consider  an  atom  that  can  diffuse  either  substitutional ly 
(via  vacancies)  or  interstitial ly.  The  flux  of  each  species  depends  on 


Implicit  in  these  arguments  is  the  assumption  that  the  ionization 
energy  stays  well  defined  even  at  high  acceptor  concentrations. 

We  present  evidence  in  Section  D  that  this  is  indeed  the  case  at 
diffusion  temperatures. 
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EFFECTIVE  DIFFUSION  COEFFICIENT  CM^SEC" 


o 


I 


Fig.  10  Zinc  diffusion  in  GaAs  at  lOOO^C  where  the  solid  line  demonstrates 
the  poor  fit  with  Allen *s  theory 


-27- 


the  concentration  gradient  of  each,  and  an  effective  diffusion  co¬ 
efficient  can  be  calculated  using  Eq.  (3)  or  (4).  In  the  following 
we  will  use  Eq.  (4)  as  an  approximation  even  where  the  fraction  in  a 
given  state  is  a  function  of  concentration.  Doing  this  will  introduce 
little  error  and  will  greatly  simplify  the  considerations.  Rewriting 
Eq.  (4)  in  terms  of  the  interstitial  and  substitution  concentrations, 
and  and  diffusion  coefficients,  D^  and  D^,  we  have: 


*1  N,  +  °S 


Nj  +  Ns 


(14) 


where  the  first  term  will  dominate  if  an  interstitial  mechanism  is 

dominant  and  the  latter  will  be  important  if  a  substitutional  mechanism 

dominates.  In  most  systems  of  interest,  the  fraction  interstitial, 

“2  -8 

N^/  +  N^,  is  small  (i.e.,  10  -  10  ),  but  is  so  much  greater 

than  D^,  that  the  interstitial  mode  may  still  dominate  D.  This  is  so 
for  the  diffusion  of  copper  and  gold  in  silicon  and  germanium.  On  the 
other  hand,  for  groups  III  and  V  atoms  in  silicon  and  germanium,  the 
evidence  suggests  that  the  substitutional  mode  dominates,  which  implies 
negligible  interstitial  concentration  for  these  impurities.  In  the 
following  sections  we  develop  defect  equilibria  that  will  enable  us  to 
calculate  the  expected  concentration  dependence  of  D  when  either  of  the 
two  modes  is  dominant. 


a.  Dominant  Interstitial  Mode 


Consider  an  impurity  that  is  an  acceptor  when 
in  a  substitutional  position,  but  is  a  donor  (positively  charged)  when 
intersti tial ly  placed.  There  will  be  an  equilibrium  reaction  in  the 
crystal  between  these  two  forms,  namely 


A  substitutional  donor  can  be  treated  in  a  similar  manner. 
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and  the  Law  of  Mass  Action  will  relate  the  concentration  of  these  species 


as : 


2 

P  ^ 


where  and  are  the  substitional  and  interstitial  concentrations, 

and  p  is  the  hole  concentration.  We  could  have  included  a  vacancy 

in  the  reaction  above,  since  it  is  certainly  involved  in  the  transfer 

of  an  atom  from  one  mode  to  another,  but  this  would  merely  introduce  a 

constant  vacancy  concentration  term  into  the  equilibrium  constant  Kj* 

The  important  point  is  that  the  fraction  interstitial  is  proportional 
2 

to  P  ,  It  can  also  be  shown  in  the  same  manner  that  the  fraction  in- 

3 

terstitial  is  proportional  to  p  if  the  interstitial  atom  becomes  a 

0 

doubly  ionized  donor.  Longini  recently  made  this  suggestion  to  explain 
our  results  for  zinc  diffusion  in  GaAs. 


We  have  not  treated  the  case  of  a  diffusing^- 
ionized  interstitial  because  we  feel  that  this  interstitial  form  is 
likely  to  diffuse  much  more  slowly  through  the  lattice  and  hence,  is 
not  likely  to  dominate  the  diffusion  process.  However,  if  such  a 
species  were  important,  the  effective  diffusion  coefficient  would  be 
proportional  to  the  first  power  of  the  hole  concentration. 

The  companion  problem  of  substitutional  donor  can 
be  treated  in  the  same  way.  However,  if  the  interstitial  is  also  a 
singly  ionized  donor,  then  the  fraction  interstitial  will  be  independ¬ 
ent  of  electron  concentration.  If  the  interstitial  is  a  doubly  ionized 
donor,  D  wi  1 1  be  depressed  at  high  substitutional  donor  concentrations. 

If  the  interstitial  form  can  diffuse  significantly  in  the  un-ionized 
state  or  exist  in  acceptor  form,  then  D  wi 1 1  be  proportional  to  the  elec* 
tron  concentration  to  some  power.  The  latter  possibility  is  believed  to 
be  extremely  remote  because  of  the  large  ionic  radius  expected  of  a 
negatively  charged  interstitial. 


In  summary,  we  feel  that  the  concentration  depend¬ 
ence  of  wi 1 1  help  support  an  interstitial  process  in  the  following 
cases: 


(1)  The  effective  diffusion  coefficient  of  a  substitutional 
acceptor  square  or  cube  of  the  acceptor  (hole)  concentration. 

(2)  The  effective  diffusion  coefficient  of  a  donor  being  de¬ 
pressed  with  high  donor  (electron)  concentrations, 

b.  Dominant  Substitutional  Mode 

Now  consider  the  case  where  the  second  term  in  Eq. 
(l4)  is  important.  The  fraction  substitutional  will  be  essentially 
unity  at  all  concentrations  and  the  diffusion  coefficient  will  be 
a  function  of  impurity  concentration.  This  occurs  because  is  a 
function  of  the  vacancy  concentration,  which  in  turn  is  affected  by 
the  impurity  concentration.  We  will  assume  that  is  directly  pro¬ 
portional  to  the  vacancy  concentration.  If  the  vacancy  exists  in 
several  states  of  ionization,  the  substitutional  diffusion  coefficient 
of  the  impurity  may  consist  of  several  modes,  each  term  of  which  will 
be  proportional  to  the  concentration  of  each  vacancy  type.  Although 
we  will  retain  this  fact,  we  will  discuss  next  the  simplest  situation 
in  which  (and  hence  D)  is  a  function  of  one  type  of  vacancy  only. 


Consider  the  case  of  a  vacancy  that  ionizes  to 
become  pos i t i ve 1 y  cha rged ,  name 1 y : 


+ 

e  + 


r 


Using  the  Law  of  Mass  Action  again,  we  have: 


+ 


P 


(16) 
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where  and  N®  are  the  positively  charged  and  neutral  vacancy  con¬ 
centrations,  p  tha  hole  concentration,  and  an  equilibrium  constant 
that  is  a  function  of  the  ionization  energy  of  the  vacancy.  The 
neutral  vacancy  concentration,  nJ;  is  not  a  function  of  impurity  con- 

''  9 

centration,  at  least  to  a  first  approximation.  Hencc^  the  positively 
charged  vacancy  concentration  will  be  proportional  to  the  first  power 
of  the  hole  concentration.  If  the  vacancy  can  become  doubly  or  triply 
ionized,  the  concentration  of  these  multiply  ionized  vacanies  will  be 
proportional  to  the  second  and  third  powers  of  the  hole  concentrations. 

Then  summarizing  the  case  for  a  dominant  substi¬ 
tution  mode  diffusion  of  an  acceptor,  one  should  expect  the  diffusion 
coefficient  to  be  enhanced  to  some  power  of  the  hole  concentration  if 
positive  vacancies  are  important.  If  negatively  charged  vacancies  are 
the  only  ionized  vacancy  allowed,  the  D  should  be  depressed  with  in¬ 
creasing  acceptor  (hole)  concentration.  The  depression  of  D  wi 1 1  con¬ 
tinue  until  the  neutral  vacancies  begin  to  dominate,  when  D  wi 1 1  become 

* 

independent  of  concentration. 


The  above  arguments  will  apply  to  paired  vacancies 
(di vacancies)  as  well  as  to  simple  vacancies.  If  divacancies  are  the  pre¬ 
dominant  mode  of  diffusion  in  III-V  compounds,  the  divacancies  could 
be  supplied  by  a  reaction  such  as 


'Ga 


+  V 


As 


where  W®  is  the  symbol  chosen  for  neutral  divacancy.  The  Mass  Action 
Law  expression  for  this  reaction  is 


Exactly  complementary  arguments  apply  to  the  effect  of  high  donor 
concentrations  on  the  concentration  of  negative  and  positive  va¬ 
cancies,  which  will  be  enhanced  and  depressed/  respectively. 
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'vGa  %As  =  Vw 


(16a) 


Since  the  product  of  the  simple  vacancy  concentration  is  usually  taken 
as  constant,'^  the  equilibrium  neutral  divacancy  concentration  will 
also  be  constant.  This  implies  that  the  diffusion  coefficient  of  an 
impurity  that  Is  dominated  by  divacancy  diffusion  will  not  be  a  function 
of  the  pressure  of  either  component  over  the  compound.  The  concentra- 
tion  of  ionized  dlvacancies  will  still  be  a  function  of  hole  and  elec¬ 
tron  concentration,  of  course. 

C.  Zinc  Diffusion  in  GaAs 

The  pages  immediately  following  contain  a  discussion  which 
shows  the  results  of  diffusing  zinc  into  gallium  arsenide  with  parti¬ 
cular  emphasis  on  the  concentration  dependence  of  D.  A  critical  analysis 
of  these  data  is  also  included. 


1 ,  Experimental  Procedure 

The  GaAs  used  was  n-type  with  about  10^^  donors  cm 
Sample  sizes  ranged  from  rectangular  samples,  6  mm  x  8  mm,  to  square 
ones,  10  mm  X  10  mm.  Sample  thicknesses  before  diffusion  ranged  from 
0.5  mm  to  1  mm.  Most  slices  were  prepared  for  diffusion  with  l400 
carborundum,  and  then  fine  lapped  on  a  flat  glass  plate  with  American 
Optical  polishing  compound  309W.  This  treatment  produced  a  surface 
that  was  monitored  for  flatness  using  simple  interference  effects  by 
utilizing  a  precision  optical  flat  and  oblique  viewing.  The  surfaces 
were  flat  to  better  than  1  micron.  The  samples  were  cleaned  by  boil¬ 
ing  in  an  HCl-methanol  solution,  then  rinsed  in  deionized  water. 

All  the  diffusion  runs  were  made  in  sealed  quartz  am¬ 
poules  with  sufficient  zinc  or  zinc-gallium  alloy,  not  in  contact  with 
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the  GaAs,  to  n^aintain  equilibrium  conditions  during  the  run.  The  zinc- 
gallium  alloy  was  used  to  control  the  surface  concentration.  After  the 
sample  was  diffused  with  radioactive  Zn^^,  the  edges  were  lapped  off. 

Zinc  distribution  was  measured  by  lapping  off  successive  layers  0,5  - 
5m.  thick.  The  amount  removed  was  monitored  by  weighing  the  sample 
with  a  Mettler  microbalance.  All  lapping  was  done  by  hand  with  AO 
309W  polishing  compound  on  2-inch  dium  optically  flat  glass  or  pressed 
alumina  discs.  Occasional  checks  with  a  precision  optical  flat  were 
made  to  ascertain  that  reasonable  flatness  (better  than  1  micron)  was 
maintained  on  both  optical  flats  and  the  sample.  The  samples  were  ro¬ 
tated  often  during  lapping  to  favor  plane  parallel  material  removal. 
Continuous  monitoring  of  the  conductivity  type  of  the  sample  using  a 
thermoelectric  probe  served  as  a  check  on  the  wedging.  For  most  samples, 
the  type  change  occurred  simultaneously  across  the  whole  slice.  After 
lapping,  the  sample  was  wiped  with  a  small  piece  of  tissue  paper  placed 
on  the  optical  flat.  The  optical  flat  with  lapping  compound  and  tissue 
paper  was  then  wrapped  with  a  single  layer  of  Saran  Wrap  and  approxi¬ 
mately  centered  on  a  well  counter.  The  exact  position  of  the  radio¬ 
active  material  on  the  2-inch  plate  was  not  critical.  After  a  count 
was  made,  the  optical  flat  was  washed  in  hot  aqua  regia. 

The  specific  activity  of  the  zinc  was  about  100  mC/gm. 

This  relatively  low  activity  limited  the  concentration  determinations 

17  -3 

to  greater  than  about  5  x  10  cm  .  However,  after  this  limit  was 
reached,  lapping  was  continued  in  small  increments  until  the  p-n 
junction  was  reached.  This  junction  depth  provided  an  accurate  zinc 
concentration  determination  at  about  10^^  cm"\  The  remainder  of  the 
sample  was  n-type  in  almost  all  cases,  with  the  bulk  resistivity  only 
slightly  higher  after  diffusion.  Little  or  no  zinc  penetration  to  the 
interior  of  the  sample  was  indicated. 
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2. 


Results 


The  results  of  most  of  the  radioactive  tracer  runs  we 
have  made  are  shown  in  Table  IV  and  in  Figs,  11-14.  The  individual 
data-points  on  the  diffusion  profiles  at  the  higher  concentrations 

19 

should  be  within  5%  of  the  true  concentration.  In  the  range  of  10 
-3  l8  "S 

cm’^  and  10  cm  ^  the  uncertainties  are  estimated  to  be  about  10%  and 
50%,  respectively. 


All  the  profiles  measured  deviated  sharply  from  those 
expected  for  a  constant  diffusion  coefficient.  The  profiles  were  ana¬ 
lyzed  for  the  effective  diffusion  coefficient  of  zinc  as  a  function  of 
zinc  concentration  using  the  Boltzmann-Matano  method  discussed  by  Crank, 
at  a  given  concentration,  C^,  on  a  diffusion  profile  is  given  by; 


7 


D  =  1 _ L-.-.- 

1  2t  dC/dx  ^ 


X  dC 


Area  under  C 
2t  dC/dx 


(17) 


where  t  is  time.  The  concentration  gradient  at  and  the  area  under 
the  curve  are  easily  obtained  from  a  plot  of  C  vs  x  on  a  linear  scale. 

The  area  may  be  measured  with  sufficient  accuracy  by  counting  squares  or 
by  using  a  planimeter. 

In  principle,  the  complete  concentration  dependence  of  0 
can  be  obtained  from  a  single  profile.  However,  as  Crank  points  out, 
the  values  at  the  extremes  of  concentration  may  be  in  considerable  error, 
especially  if  the  D  variation  is  large.  The  region  near  the  surface  may 
be  in  error  because  of  out-diffusion  during  the  cooling  cycle  or  other 
effects  associated  with  the  surface.  Errors  may  be  introduced  in  the 
regions  of  low  concentration  from  the  difficulty  in  recording  the  exact 
slope  in  this  vicinity.  Diffusion  profiles  for  the  zinc  in  GaAs  system 
have  quite  high  concentration  gradients  in  the  region  near  the  p-n  Junc¬ 
tion  and  the  D  values  measured  are  invariably  higher  than  expected.  This 
may  be  from  mild  wedging  of  the  samples  during  analysis,  but  could  also 
be  from  a  non-equilibrium  effect  in  this  vicinity. 
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Table  IV 


Z i nc  Diffusion  in  GaAs 


Run 

Percent  Zn 
in  Source 

Temperature 

rc) _ 

1 

100.0 

o 

O 

o 

a\ 

2^ 

100.0 

900* 

3 

1.0 

VO 

o 

o 

o 

4 

* 

0.25 

900* 

5 

3.0 

900* 

6 

100.0 

750" 

7 

5.0 

750* 

8 

2.2 

750° 

9 

1.0 

750° 

10^^ 

100.0 

692° 

2.0 

1.5 
87.5 
166.0 

16.0 

5.5 
25.0 
25.0 

303.0 

8.0 


p-n  Junction 


208 

71 

60 

23 


75 

13-15 

9-10 

15-16 

14 


Amount  of  zinc  limited  to  provide  130  torr  vapor  pressure;  incremental 
resistivity  also  measured  to  determine  electrically  active  concentration. 
'it 

Iso-concentration  run  diffusing  radio-zinc  into  already  zinc-doped  sample* 
Both  sides  of  sample  analyzed  to  check  for  orientation  effects. 
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Fig.  11 


Fig.  12  Iso- concentration  run.  Surface  concentration  with  tracer 
zinc  maintained  at  same  concentration  as  non-tracer  bulk 
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Fig.  14  Zinc  diffusion  profile  in  gallium  arsenide  at  692*C 
showing  absence  of  orientation  effects 


For  the  reasons  discussed  above,  we  have  ignored  values 
of  0  taken  from  regions  very  near  the  surface  and  near  the  Junction. 
Results  of  the  analysis  of  the  profiles  at  900*C  are  shown  in  Fig.  15. 
Also  shown  is  our  analysis  of  the  Gunnel  1  and  Gooch^  profiles  at 
this  temperature  and  the  results  obtained  at  Stanford  by  Pearson  and 
co*-workers.  Agreement  among  these  three  laboratories  is  striking 
for  runs  made  under  similar  conditions.  For  example,  note  the  run 

made  at  the  highest  concentration  by  each  group,  in  which  an  apparent 

—  19  *3 

maximum  in  D  is  recorded  by  each  group  at  about  7  x  10  cm  •  The 

effective  diffusion  coefficients  at  a  given  concentration  agree  to 

within  about  10%,  which  is  extremely  good  for  such  a  complex  system. 

This  agreement  in  data  is  especially  surprising  in  view  of  the  quite 

different  experimental  technique  employed.  They  used  a  two<-temperature 

method  to  control  their  surface  concentration  instead  of  the  zinc- 

gallium  source.  For  runs  made  at  lower  concentrations,  the  same  sort 

of  effect  was  observed,  but  the  apparent  maxima  occurred  at  a  lower 

concentration. 


These  efforts  indicate  that  the  measured  effective  dif¬ 
fusion  coefficient  is  not  only  a  function  of  concentration,  but  also  a 
function  of  some  other  variables,  such  as  time  and  distance.  The 
method  used  in  the  analysis  of  D  is  good  only  if  D  is  a  function  of 
concentration  only.  Since  this  is  obviously  not  the  case  for  Zn  in 
GaAs,  we  have  developed  a  technique  which  subordinates  these  diffi¬ 
culties.  It  involves  diffusing  with  radioactive  zinc  into  material 
already  doped  with  non- radioactive  zinc.  Diffusion  conditions  are 
set  up  so  the  expected  surface  concentration  is  identical  to  the  zinc 
concentration  already  present  in  the  bulk.  We  call  this  an  iso¬ 
concentration  diffusion,  since  the  indiffusing  radioactive  zinc  always 
finds  itself  in  a  constant  total  zinc  concentration  ambient.  Because 
there  are  no  complicating  zinc  concentration  gradients,  the  expected 
profile  is  of  a  system  with  constant  diffusion  coefficient  and  constant 
surface  concentration,  namely  a  complementary  error  function. 


EFFECTIVE  DIFFUSION  COEFFICIENT  CM^  SEC 


Fig.  15  Zinc  diffusion  in  GaAs  at  900°C 
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Figure  12  shows  the  results  obtained  by  diffusing  radio¬ 
active  zinc  from  an  alloy  of  Zn^^  and  gallium  (3%  Zn) ,  known  from  pre- 

19  -3 

vious  data  to  control  the  surface  concentration  at  about  2  x  10  cm  , 
into  GaAs  already  doped  during  crystal  growth  at  this  same  level  with 
non- radioactive  zinc.  The  expected  erfc  distribution  is  followed  quite 
closelyt  and  the  diffusion  coefficient  is  about  a  factor  of  2  lower  than 
that  determined  by  the  previous  method. 

An  alternative  method  for  obtaining  a  non-radioactive 
sample  at  the  desired  concentration  is  by  a  lengthy  pre-diffusion  pro¬ 
cess  under  exactly  the  same  diffusion  conditions  as  the  proposed  tracer 

run.  This  method  has  been  successfully  applied  by  the  Stanford  group 

—  20  -3 

and  a  0  value  obtained  at  a  zinc  concentration  of  3  x  10  cm  ,  of 

—8  2  10 

about  6  X  10  cm  /sec  .  This  result  is  much  higher  than  the  values 
obtained  by  analyzing  the  normal  diffusion  profiles,  so  the  apparent 
depression  of  D  at  extremely  high  concentrations  probably  results  from 
kinetic  effect  such  as  non-equilibrium  vacancy  and/or  interstitial 
conditions.  Utilizing  these  facts,  we  have  drawn  a  curve  to  indicate 
where  we  think  the  ‘^true"  effective  diffusion  coefficient  should  fall. 
This  curve  should  be  regarded  as  quite  tentative,  of  course,  since  it 
heavily  favors  only  two  points,  namely  the  iso-concentration  runs. 

Also  shown  in  Fig.  11,  on  run  2,  are  the  results  of  an 
incremental  resistivity  method^ ^  which  was  undertaken  in  conjunction 
with  radioactive  sectioning.  This  was  done  to  determine  the  extent 
to  which  the  zinc  was  electrically  active.  To  convert  from  resis¬ 
tivity  to  a  carrier  concentration,  we  used  the  combined  data  of  Haisty 
12 

and  Kellett  shown  in  Fig.  16.  The  agreement  is  within  about  15%  at 
most  concentrations,  which  is  within  the  probable  error  of  the  resis¬ 
tivity  measurements.  Zinc  is  apparently  completely  electrically 
active  at  room  temperature  in  these  concentrations.  At  still  higher 

concentrations  evidence  shows  this  is  not  so.  For  example,  with  zinc 

20  -3 

concentration  of  about  3  x  10  cm  ,  the  hole  concentration  in  one 

20  -3 

sample  was  only  about  1  x  10  cm  .  We  take  this  as  evidence  that 
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at  these  very  high  concentrations  the  zinc  precipates  significantly 
during  the  cooling  cycle. 

To  prove  that  precipitation  could  actually  occur,  even 
at  the  relatively  low  temperature  of  220^0,  we  continuously  monitored 
the  resistivity  of  one  sample  for  several  months  at  this  temperature. 
Contacts  were  applied  using  a  high  temperature  solder.  These  results 
are  shown  in  Fig.  17^  with  the  average  hole  concentration  being  esti¬ 
mated  from  the  resistivity  at  220^C. 

-n  -14  2 

A  diffusion  coefficient  of  10  ^  to  10  cm  /sec  can  be 

estimated  at  220^C  for  the  initial  part  of  the  precipitation  by  making 

certain  simplifying  assumptions  concerning  the  shape  of  the  precipitate 

and  the  concentration  of  precipitation  centers.  The  process  slowed 

somewhat  as  the  average  active  zinc  concentration  decreased,  but  not 

as  much  as  expected.  Several  zinc-doped  crystals  in  the  range  10^^ 

“3  1 8  —3 

cm  ^  to  10  cm  ^  were  monitored  in  the  same  manner  at  220°C  and  only 
negligible  increases  in  the  resistivity  were  observed  (less  than  3% 
in  64  days). 


The  implications  of  the  above  measurements  are  clear. 

As  an  emitter  in  a  transistor  or  in  tunnel  diodes,  some  precipitation 

19  -3 

of  zinc  should  be  expected  at  concentrations  above  10  cm  .  At 
l8  *"3 

concentrations  of  10  cm  ^  and  lower,  the  zinc  should  remain  in 
solid  solution  for  extended  periods,  even  at  reasonably  high  opera¬ 
ting  temperatures. 

Results  at  7S0^C  shown  in  Figs.  13  and  18  are  similar  to 

those  at  900^C.  The  apparent  maximum  in  D  occurs  at  about  4  x  10^^ 

-3  l4 

cm  .  An  additional  diffusion  coefficient  has  been  provided  by  O^Asaro 

at  an  average  zinc  concentration  of  about  7  x  10^^  cm*^.  At  this  con- 

-14  2  — 

centration  he  obtained  a  value  of  about  3  x  10  cm  /sec  for  0,  off 
scale  on  Fig.  18^  but  consistent  with  the  dotted  line  representing  our 
estimate  of  the  effective  diffusion  coefficient. 
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Fig.  17  Effect  of  heat  treatment  at  220®C  on  heavily  zinc  doped  GaAs  demonstrating  precipitation. 


piffosion  coefficient  cm^ 


The  results  at  lOOO'^C  were  considered  in  the  theory 
section,  where  it  was  shown  that  Allen's  suggestion  that  the  variation 
in  D  was  due  to  the  different  diffusion  rates  of  ionized  and  un-ionized 
zinc  was  not  valid.  Most  the  data  shown  is  from  reference  1^  al¬ 
though  one  set  of  data  at  high  concentration  1$  from  the  Stanford 
group At  low  concentrations^  the  recent  data^^  are  shown.  The 
dotted  line  again  is  a  tentative  estimate  of  the  effective  diffusion 
coefficient. 


Our  analysis  of  the  Gunnell  and  Gooch  data  taken  with 
excess  arsenic  in  the  ampoule  at  lOOO^C  Is  also  shown  in  Fig,  10.  The 
effective  diffusion  coefficient  is  apparently  markedly  reduced  under 
these  conditions  near  the  surface,  but  not  affected  significantly 
deeper  in  the  sample.  We  are  not  sure  if  this  is  an  equilibrium  or 
a  kinetic  effect. 

An  apparent  activation  energy  for  diffusion, for 

fl 

zinc  in  GaAs  can  be  obtained  from  the  three  sets  of  data  £t  ^  given 
concentration.  The  apparent  activation  energy  is  significantly  reduced 
at  high  concentrations.  The  expression  relating  D  to  temperature  is 


D  »  D  exp 

oa  ^ 

[-V"]  • 

(18) 

where  both  D  and  Q.  are  a  function  of  concentration.  The  values  ob¬ 
oe  a 

tained  for  Q  and  D  at  several 
a  oa 

concentrations 

fol low: 

Zinc  Concentration 

Q,  (eV) 

2 

DqA  (cm  /sec) 

3  X  1o'^ 

0.6  ±  0.2 

6  X  10“^ 

1  X 

1.0  t  0.2 

2  X  10‘® 

1  X  lo’® 

1.6  1  0.2 

4  X  10’® 

Intrinsic 

2.3  ±  0.2 

5  X  10"^ 
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Values  quoted  for  intrinsic  GaAs  were  obtained  by  extrapolation  of  the 
experimental  data  to  zinc  concentrations  below  the  intrinsic  electron 
concentration  at  a  given  temperature  where  D  is  expected  to  become  in¬ 
dependent  of  zinc  concentration.  It  should  be  emphasized  that  the 
estimated  0  at  220^C  obtained  from  the  precipitation  experiment  leads 
to  an  activation  energy  of  about  0.7  ^  0.2  at  high  zinc  concentrations. 

This  is  consistent  with  the  Q  estimated  above  in  this  range  of  concen- 

a 

trations. 


A  few  of  the  tracer  runs  made  were  analyzed  on  both  sides 
of  the  sample.  A  run  at  692^C>  shown  in  Fig.  14,  demonstrating  that  there 
is  no  significant  effect  of  crystal  orientation  on  D.  Although  none 
was  expected,  there  was  a  possibility  that  the  surface  equilibrium  con¬ 
ditions  would  be  different  on  opposite  sides  of  a  <11 1>  slice,  since  a 
different  type  atom  predominates  on  each  side.  No  such  effect  was  ob¬ 
served  in  this  run,  however. 

As  an  engineering  aid,  we  have  compiled  most  of  the 

existing  zinc  data  as  to  p-n  Junction  depth  and  have  shown  them  in  Fig. 

19  plotted  against  zinc  surface  concentration.  Surface  concentration 

can  be  predicted  accurately  from  the  diffusion  conditions  from  data  to 

be  presented  in  Section  0.  The  junction  depth  is  that  expected  for  an 

17-3  4 

n-type  bulk  concentration  of  10  cm  and  for  a  time  of  10  seconds. 

The  values  were  calculated  by  assuming  the  Junction  depth  advanced  as 
the  square  root  of  time.  This  is  a  valid  assumption  even  for  a  sharp 
concentration  diffusion  coefficient,  but  assumes  that  0  is  a  function 
of  concentration  only.  The  anomalous  effects  discussed  earlier  in¬ 
dicate  that  D  is  evidently  a  function  of  time  and  distance  in  this 
system,  but  the  data  of  Fig.  19  can  be  used  as  a  first  approximation. 

The  inflection  points  in  both  the  900*^0  and  lOOO^C  are  believed  due 
to  the  kinetic  effects  discussed  earlier.  In  most  runs  reported,  it 
is  likely  that  the  surface  concentration  had  reached  a  constant  value. 
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1000 


ZINC  SURFACE  CONCENTRATION 


Fig.  19  P-N  junction  depth  for  zinc  diffusion  in  GaAs  for  a  normalized 
time  of  10^  seconds  and  10^  'cm*^  bulk  donors 


We  have  evidence^  however,  that  the  attainment  of  an  equilibrium 
surface  concentration  is  rather  sluggish,  so  short  diffusion  times 
will  generally  lead  to  shallower  Junctions  than  those  predicted  here, 

3.  Discussion  of  Results 


The  important  experimental  facts  to  be  explained  for 
the  zinc  in  GaAs  system  at  any  given  temperature  are: 

(1)  ^  is  proportional  to  where  n  ranges  from  2  to  2,5# 

(2)  An  apparent  maximum  in  D  occurs  when  the  diffusion  profiles 
are  analyzed. 

(3)  High  arsenic  pressures  tend  to  depress  D  near  the  surface. 
These  experimental  facts  can  be  explained  assuming  ei ther  a  dominant 
interstitial  diffusion  mode  or  a  dominant  substitutional  mode. 

For  an  interstitial  process:  Fact  (!)  can  be  explained 
assuming  that  the  zinc  interstitial  is  a  singly  (sometimes  doubly) 
ionized  donor.  Fact  (2)  can  be  justified  on  the  assumption  that  the 
high  diffusion  coefficients  in  the  interior  are  caused  by  the  higher 
than  equilibrium  "fraction  interstitial"  (which  dominates  D,  as  shown 
in  the  theory  section),  since  the  vacancies  in  this  region  must  be 
generated  at  dislocation,  a  sluggish  process.  The  region  near  the 
surface  is  nearer  the  equilibrium  fraction  interstitial,  since  the 
surface  is  a  good  source  of  incoming  vacancies,  so  D  is  lower  near 
the  surface  than  in  the  interior.  Fact  (3)  can  be  explained  by  assum¬ 
ing  the  gallium  vacancies  are  enhanced  in  the  presence  of  high  arsenic 
pressure.  If  the  interstitial  zinc  tends  to  react  predominantly  with 
gallium  vacancies,  then  the  enhanced  gallium  vacancy  concentration  would 
causa  the  fraction  Interstitial  (and  hence  D)  near  the  surface  to  be 
depressed.  In  the  interior  of  the  sample,  D  will  remain  high  because  of 
the  dislocation  generated  arsenic  and  gallium  vacancies. 
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For  a  substitutional  process  the  arguments  are  the 
following.  Fact  (1)  requires  that  the  vacancies  be  doubly  (or  sometimes 
perhaps  triply)  ionized  donors.  These  must  be  divacancies,  since  simple 
gallium  vacancies  are  inconsistent  with  Fact  (3),  as  will  be  shown. 

Fact  (2)  can  be  explained  by  postulating  a  flux  of  dislocation  generated 
vacancies  from  the  interior  into  the  heavily  p-type  region.**"  This  will 
cause  higher  D  values  than  typical  of  equilibrium  at  depths  near  the 
junction  since  a  flux  of  vacancies  outward  causes  a  net  flux  of  atoms  inward. 
The  vacancy  concentration  near  the  surface  may  be  depressed  by  the  incoming 
zinc^  because  the  zinc  can  enter  the  crystal  only  at  a  vacant  site. 

This  would  explain  the  low  D  values  near  the  surface.  At  long  times 
the  net  flux  of  zinc  atoms  into  the  crystal  would  finally  decrease  and 
incoming  vacancies  could  maintain  their  equilibrium.  The  sluggish 
approach  to  the  equilibrium  surface  concentration  mentioned  earlier  may 
also  result  from  this  effect. 

Fact  (3)  cannot  be  explained  under  the  assumption  that 
simple  gallium  vacancies  are  the  main  mode  of  diffusion  for  zinc.  This 
is  evident  because  excess  arsenic  should  enhance  the  gallium  vacancy 
concentration,  which  would  in  turn  enhance  D.  On  the  other  hand,  the 
way  in  which  the  di vacancy  concentration  could  be  reduced  is  a  non* 
equi I ibrium  effect  due  to  a  sluggish  gallium  vacancy  equilibrium. 

The  arsenic  vacancy  concentration  will  be  sharply  reduced  in  an  excess 
arsenic  atmosphere.  The  gal  1 ium  concentration  should  increase  propor* 
tionally  to  keep  the  gal  1 ium*arsenic  vacancy  product  constant.  If  the 
latter  process  is  sluggish,  the  divacancy  reaction  will  not  be  supplied 


An  alternative  explanation  for  Fact  (1)  involves  a  singly  charged 
di vacancy  whose  ionization  energy  decreases  with  increasing  acceptor 
concentration.  This  could  be  expected  in  view  of  the  large  orbital 
of  the  electron  in  a  material  like  GaAs.  The  electronic  orbit  of  a 
donor  would  contain  several  acceptors  at  these  high  zinc  concentrations. 

^  Some  very  unusual  data  are  presented  in  a  later  section  about  cadmium 
diffusion  that  seem  to  substantiate  the  outward  vacancy  flux  concept. 
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properly  so  the  concentration  of  the  latter  will  decrease,  as  will  0« 
Near  the  Junction,  D  will  still  be  high  because  of  bulk  generated 
vacancies , 


In  looking  at  the  two  choices  available,  one  is  tempted 
to  accept  the  interstitial  argument  and  reject  the  substitional  because 
of  the  greater  simplicity  of  the  arguments  in  the  interstitial  case. 
However^  in  addition  to  the  cadmium  data  which  substantiates  the 
vacancy  flux  proposal^  there  are  three  other  semiconductor  systems 
that  can  be  explained  very  nicely  with  substitutional  models^  whereas 
the  interstitial  model  is  inconsistent  in  almost  every  case.  In  silicon 
the  evidence  is  that  divacancies  do  indeed  dominate  diffusion  at  high 
temperatures^  and  furthermore^  that  they  can  be  singly  ionized  donors > 
neutral^  or  singly  or  doubly  ionized  acceptors. 

4.  Summary 

We  believe  that  the  results  for  zinc  diffusion  in  GaAs 
are  ambiguous  for  determining  a  diffusion  mechanism  from  the  observed 
concentration  dependence  of  the  diffusion  coefficient.  The  preponder¬ 
ance  of  data  in  other  semiconductor  systems  seems  to  favor  a  substitu¬ 
tional  mechanism  for  them,  but  this  does  not  shed  any  direct  light  on 
the  GaAs  problem.  We  feel  that  self-diffusion  data  for  gallium  and 
arsenic  in  heavily  doped  GaAs  will  perhaps  resolve  the  question. 

D.  Zinc  Surface  Concentration  Control 


1,  Vapor-Solid  Equilibria 

The  surface  concentration  of  an  impurity  in  a  semicon¬ 
ductor  during  diffusion  can  be  specified  by  controlling  the  vapor 
pressure  of  the  impurity  over  the  semiconductor.  We  will  develop 
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equations  relating  the  equilibrium  surface  concentration  to  the  vapor 
pressure. 


Consider  an  acceptor  chat  exists  monatomically  in  the 
vapor  that  is  in  equi librium  with  a  semiconductor.  The  reaction  involv¬ 
ing  transfer  of  the  atom  from  the  gas  to  the  solid  is 


» 


where  the  atom  in  the  semiconductor  is  assumed  to  be  un-ionized  until  a 
subsequent  reaction  occurs,  namely: 

"a 

A®  ^  a”  + 


The  Mass  Action  Law  applied  to  the  first  reaction  gives 

"a  -Va  •  ('« 

where  is  the  un-ionized  acceptor  concentration  in  the  solid,  P. 
the  partial  pressure  of  the  acceptor  and  an  equilibrium  constant. 
The  second  reaction  above  has  already  been  discussed  and  the  resulting 
mass-action  expression  given  in  £q.  (5)  as 


(5) 


where  p  is  the  hole  concentration  and  can  be  calculated  using  (6). 


By  utilizing  these  expressions  and  the  charge  neutrality 
and  conservation  of  species  as  before,  the  following  expression  can  be 
developed  (assuming  no  compensating  donors  are  present): 


Wf< 


(20) 
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where  is  the  equilibrium  acceptor  concentration  in  the  solid,  and 
nj  the  intrinsic  electron  concentration.  The  left  and  right  terms 
are  the  un-ionized  and  ionized  acceptors,  respectively. 

At  low  partial  pressures  where  the  acceptor  concen* 
tration  in  the  solid  is  much  less  than  n^,  the  above  expression  shows 
a  linear  dependence  of  on  pressure,  namely 


At  higher  pressures  where  is  much  greater  than  n^,  Eq.  (20)  be¬ 
comes: 

^A  “  '^g’’a  *  (’^g’^A^’a) 

where  the  square  root  dependence  dominates  at  intermediate  concent ra* 
tions  and  the  linear  term  becomes  important  at  very  high  acceptor  con¬ 
centrations* 


We  have  plotted  the  zinc  surface  concentrations  of  al¬ 
most  all  the  runs  made  in  GaAs*  as  a  function  of  the  pressure  of  zlnc» 
under  the  assumption  that  Raoul t's  law  is  valid  for  the  zinc-gallium 
alloy  diffusion  source*  These  data  are  shown  in  Figs*  20  and  21»  with 
data  from  the  Stanford  group  and  from  reference  1*  Pressure  in  the 
latter  was  calculated  from  the  expected  vapor  pressure  of  elemental 
zinc  which  was  placed  at  the  cold  end  of  the  tube.  The  consistency 
of  these  points  with  our  data  at  900^C  is  especially  gratifying#  con¬ 
sidering  the  different  experimental  conditions. 


Raoult*s  law  states  that  the  partial  pressure  of  the  solute  (zinc) 
in  a  solution  is  ^A  fraction  of  the 

acceptor  and  is  the  vapor  pressure  of  the  pure  acceptor  at  the 

temperature  in  question* 
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Fig.  20  Effect  of  zinc  partial  pressure  on  surface  concentration  with  theoretical 
curves  assuming  zinc  ionization  energy  of  ,08ev 


The  theoretical  curves  shown  are  calculated  from  Eq. 

(22)  assuming  an  acceptor  ionization  energy  of  .08  eV»  which  is  a 
theoretical  value  assuming  a  simple  hydrogenic  model.  This  choice 
is  somewhat  arbitrary,  because  there  are  no  reliable  published  values. 
The  data  can  be  fitted  reasonably  well  with  any  ionization  energy  in 
this  range.  K  at  each  temperature  was  calculated  using  the  band 

e  6 

structure  parameters  of  Ehrenreich  as  discussed  in  the  Theory  section. 

The  equilibrium  constant  was  calculated  to  give  the 
best  fit  with  the  data  at  each  temperature.  These  values  of  are 
shown  in  Fig.  22,  with  two  points  calculated  from  single  runs  utiliz¬ 
ing  elemental  zinc.  The  temperature  dependence  of  is  given  by: 

Kg  =  1.5  X  exp  [.53/kT]  ,  (23) 

where  the  units  of  are  cm  ^  torr"*  and  the  heat  of  reaction  is  in 
electron  volts  (with  an  estimated  uncertainty  of  only  0.04  eV) .  By 
utilizing  this  relation  for  and  Eq.  (23)  (general)  or  (22)  (extrin¬ 
sic),  the  equilibrium  surface  concentration  expected  of  almost  any  ex¬ 
perimental  condition  should  be  predictable  to  within  about  10%  accuracy. 
This  is  based  on  data  from  Fig.  20.  Only  one  point  of  the  23  shown  is 
outside  this  error. 

Note  that  the  heat  of  transfer  for  an  atom  in  the  vapor 
going  to  an  un-ionized  atom  in  the  solid  is  negative.  This  is  evident 
in  Fig.  20  because  at  a  given  pressure  the  surface  concentration  de¬ 
creases  with  increasing  temperature. 


. — ^  *  . . 

One  torr  is  a  mm  of  Hg  of  pressure. 
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These  data  are  believed  to  be  the  first  experimental 
verification  of  the  concept  proposed  by  Reiss,  that  vapor-solid 
equilibria  are  governed  by  the  un-ionized  impurity  concentration  in 
the  solid  rather  than  by  the  total  impurity  concentration.  In  most 
other  diffusion  systems,  either  the  partial  pressure  of  the  impurity 
is  unknown  or  the  ionization  energy  is  a  function  of  concentration. 
This  latter  effect  causes  the  energy  level  to  merge  with  the  nearest 
band  at  high  concentrations  and  the  impurity  is  completely  ionized  at 
essentially  all  concentrations.  When  this  occurs^  the  theory  is  not 
fol lowed. 


Some  data  of  the  same  nature  reported  here  for  the  zinc 

in  InSb  system  is  in  poor  agreement  with  this  theory.  This  is  not 

unexpected  in  view  of  the  larger  expected  orbitals  in  this  system.  An 

order  of  magnitude  criterion  for  the  concentration  at  which  merging 

should  occur  is  given  by  calculating  the  concentration  at  which  the 

donor  or  acceptor  orbitals  begin  to  overlap.  On  a  simple  hydrogenic 

2 

model  the  radius  of  the  orbital  is  given  by  ao^  mo/irf^  where  ao  is 

the  radius  of  the  electron  orbit  of  the  hydrogen  atom  (0.53A®),  e 

the  dielectric  constant  and  m©  and  m*  are  the  electronic  and  effective 

masses,  respectively.  For  acceptors  in  InSb,  merging  of  the  energy 

level  with  the  valence  band  on  this  model  should  occur  at  about 
1 8  -3 

3  X  10  cm  so  it  is  not  surprising  that  the  theory  does  not  apply 
to  this  system,  since  all  our  data  are  for  concentrations  above  this 
value. 


The  same  calculation  for  acceptors  in  GaAs  yields  a 
20  -3 

value  of  4  X  10  cm  ,  so  the  theory  should  fit  at  concentrations  be¬ 
low  this  value.  The  good  fit  with  theory  is  good  evidence  that  the 
energy  level  of  zinc  remains  well-defined  at  these  high  tempera¬ 
tures.  We  stress  this  latter  point  because  there  is  good  evidence 
that  at  room  temperature  the  level  does  merge  with  the  valence 
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12 

band.  This  is  evident  from  the  Hal)  data  of  Haisty  and  Kellett 

shown  in  Fig.  )(>,  which  with  other  datd>  shows  zinc  to  be  completely 

20  -3 

ionized  to  at  least  10  cm  .  This  seems  to  indicate  that  the  phe¬ 
nomenon  of  impurity  level  merging  is  more  likely  to  occur  (that  is,  at 
a  lower  concentration)  at  low  temperatures.  We  feel  this  merging  is 
due  to  the  increased  degeneracy  effects  at  lower  temperatures  due  to 
the  lower  density  of  states.  An  interesting  consequence  of  this  effect 
is  that  the  hole  or  electron  concentration  of  a  heavily  doped  sample 
may  actually  decrease  as  the  temperature  is  raised  while  the  energy 
levels  *'unmerge"  at  high  temperatures. 


2,  Zinc  Solid  Solubility  and  Segregation  Coefficient 


The  segregation  coefficient,  k,  is  defined  (for  an 
acceptor,  for  example)  as 

N. 

k - ^  .  (24) 

where  is  the  acceptor  concentration  in  the  solid  and  the  ac¬ 
ceptor  concentration  in  the  liquid.  We  can  introduce  a  liquid  into 
our  vapor-solid  equilibria  if  we  imagine  that  the  surface  of  the  semi¬ 
conductor  contains  a  liquid  layer  that  is  in  equi librium  with  both  the 
acceptors  in  the  vapor  and  the  solid.  The  pressure  of  the  acceptor  is 
related  to  the  acceptor  concentration  in  the  liquid  by 


P 


A 


A 


(25) 


where  is  the  pressure  over  pure  acceptor,  the  activity  coeffi¬ 
cient  for  the  acceptor  in  the  liquid,  and  the  mole  fraction  of  the 

L  L  ^ 

acceptor  in  the  liquid,  (N.  ■  No  X.  where  No  is  the  total  atom  con> 

3  ”  ^22 

centration  per  cm  ,  which  is  4.4  x  10  for  GaAs). 
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The  segregation  coefficient  in  terms  of  the  liquid  con¬ 
centration  can  be  obtained  from  (20) »  (24),  and  (25)  and  the  result  is 


''g  V 


1  + 
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"l  ) 


] 


(26) 


where  the  activity  coefficient  (usually  of  order  unity)  can  be  inserted 
with  each  P^*i  if  necessary. 

In  terms  of  the  acceptor  concentration  in  the  solid,  tha 
complete  expression  is  rather  complex  and  will  not  be  given  here.  At 
acceptor  concentration  below  n^,  namely  intrinsic  conditions,  the  ex> 
press ion  reduces  to 

end  under  extrinsic  conditions  (N^  »  n^)  the  expression  for  k  is: 


The  ecceptor  concentration  must  vanish  at  the  melting 


point,  so  we  are  in  the  intrinsic  range  near  the  melting  point,  and 
Eq.  (27)  applies.  To  calculate  the  segregation  coefficient  near  the 
melting  point,  all  we  need  is  a  value  of  Kq  and  P^*,  since  and  n^ 
are  both  calculable,  A  value  of  Kq  can  be  obtained  et  the  melting 
point  of  GaAs  (1240*C)  from  Eq.  (23),  and  the  pressure  of  pure  line 
at  this  temperature  is  available  from  standerd  vapor  pressure  data. 

The  value  of  nj  is  1,5  x  lo'®  cm"^,  P^*  is  9  x  10®  torr,  Kq  is 

8.7  X  lo'®  cm“®  torr"',  and  is  4.6  x  lo'^  cm"®  for  an  ioniiation 
energy  of  .08  eV.  Putting  these  values  in  (27),  we  obtain  a  segregation 
coefficient  of  0.55  at  the  melting  point.  The  only  assumption  we  have 
made  is  that  the  liquid  behaves  ideally  (v^  *  1).  Tha  only  published 


(28) 
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value  of  k  for  zinc  in  GaAs  is  by  Whelan*  who  says  it  is  somewhere 

in  the  range  of  0.3  to  0.9.  Our  calculated  value  falls  in  the  middle 
of  this  range.  A  more  accurate  statement  of  our  results  is  that 
k  ■  O.SSy^;  so  knowledge  of  would  permit  a  more  accurate  value  of 
k  to  be  determined. 


The  solid  solubility  of  zinc  in  GaAs  can  be  calculated 

by  relating  the  concentration  in  the  solid  to  that  in  the  liquid  using 

^fc 

Eqs.  (20)  and  (25)  so  we  have 


^  -  '^G  ''a  Pa" 


1  + 


L  (''g  "a  ^a'  ’’a*  "l') 
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(29) 


The  temperature  dependence  of  all  these  terms  except  is  known.  This* 

too*  can  be  calculated  if  the  heat  of  fusion  of  the  semiconductor  in  known 

and  by  making  certain  simplifying  assumptions  concerning  the  ideality 
1 8 

of  the  liquid.  This  is  somewhat  complex  unless  k  «  1,  so  we  will 
show  the  data  which  we  think  best  represent  the  solid  solubility  and 
will  indicate  the  expected  temperature  dependence  as  a  dotted  line. 

This  is  done  in  Fig.  23*  where  we  have  also  included  data  from  Pearson*  et 
al.^^  The  solid  solubi 1 ity  was  obtained  by  extrapolating  the  diffusion 
profiles  to  the  surface. 


E.  Cadmium  Diffusion  in  GaAs 


Diffusion  of  radioactive  Cd^^^  in  gallium  arsenide*  as  a  p*type 
impurity*  was  investigated.  A  complication  was  observed  from  the  simuN 
taneous  diffusion  of  In^^^*  which  is  a  decay  daughter  of  Cd^^^.  It  was 


The  activity  coefficient  can  be  inserted  in  this  and  subsequent 
relations  by  multiplying  into  all  terms. 
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found  that  the  indium  had  a  diffusion  coefficient  slightly  higher  than 

the  cadmium  at  IIOO^C.  Figure  24  demonstrates  these  results,  where  the 

curve  showing  a  value  for  the  diffusion  coefficient  of  3-5  x  lO"^^ 

cm^  sec"^  is  attributed  to  the  diffusion  of  the  Cd*^^,  These  two  curves 

were  obtained  by  evaluating  opposite  sides  of  a  single  wafer,  diffused 

with  Cd^^^  for  10  hours  at  1100*C,  at  two  different  times.  The  curve 

showing  the  higher  "apparent"  surface  concentration  and  faster  diffusion 

was  obtained  immediately  after  the  diffusion  run.  The  second  curve  was 

obtained  on  the  opposite  surface  two  days  later.  The  difference  was 

attributed  to  the  fact  that  after  two  days,  the  concentration  of  the 

4.5”*hour  half-life  Zn^^^  was  down  by  3  orders  of  magnitude,  so  only  the 

53-hour  Cd^^^  was  seen.  It  Is  assumed,  therefore,  that  the  second  curve 

with  a  D  ■  2  X  lO"**  cm^  sec*^  is  due  to  Cd^'^  only.  This  value  of 

diffusion  coefficient  Is  lower  by  a  factor  of  two  than  that  reported 
19 

by  Goldstein  ;  however,  his  diffusion  conditions  were  somewhat  differ¬ 
ent. 


That  the  aforementioned  conclusions  are  probably  correct  was 
later  substantiated  in  a  diffusion  run  by  gamma  ray  spectroscopy.  This 
analysis  showed  that  the  sample  contained  in  the  regions  near  the  sur¬ 
face  an  inordinate  amount  of  In^^^  that  could  not  be  attributed  to  the 
normal  equilibrium  decay  scheme  of  Cd^^^,  The  analysis  showed  one 
other  point  of  interest:  there  was  no  detectable  concentration  of 
either  In^^^  or  Cd^^^  in  the  bulk  of  the  sample  beyond  30  microns 

penetration.  This  is  in  disagreement  with  the  results  of  Gunnell  and 
20 

Gooch,  which  were  possibly  obtained  under  somewhat  different  experi¬ 
mental  procedures.  Figure  25  shows  a  concentration  profile  for  Cd^^^ 
after  diffusion  at  lOOO^C  for  53  hours.  The  diffusion  coefficient  is 

8  X  10'  cm  sec  ,  which  is  again  lower  by  a  factor  of  four  than  that 

20  19 
reported  by  Gunnell  and  Gooch  and  a  factor  of  ten  lower  than  Goldstein 

Since  our  diffusion  was  done  under  conditions  which  gave  lower  surface 

concentrations,  it  could  be  noted  that  cadmium,  like  zinc,  is  sharply 

concentration  dependent. 
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Fig.  24  10  houra  at  1100°C  with  0.2  atm.  cadmium^^5  showing  transient  effects 

,  due  to  indium^ 


F.  Manganese  Diffusion  In  GaAs 


Diffusion  of  manganese  as  an  acceptor  in  gallium  arsenide  was 
investigated.  Several  diffusion  runs  were  made  under  different  condi* 
tions  and  the  resulting  manganese  distributions  determined  using  in¬ 
cremental  sheet  resistance  techniques.  The  results  were  in  agreement 
with  those  expected  from  the  typical  complementary  error  function 
distribution^ 


C(x)  -  C  (1-erf  - - - )  .  (30) 

From  Eq.  (30)/  values  of  0  were  calculated  for  the  various  diffusion 
conditions  and  are  plotted  vs  IOOO/TC’K)  in  Fig.  26.  The  data  may  be 
represented  by 

0  -  8.5  X  10“^  c-»-7/kT 

Figure  27  shows  the  dependence  of  the  surface  concentration  of  manganese 
on  temperature.  ,  ^ 

G.  Sulfur  Diffusion  in  GaAs 


The  study  of  diffusion  of  sulfur  into  gallium  arsenide  was  greatly 
hindered  by  reaction  of  the  diffusant  vapor  with  the  sample  at  the  necessary 
high  temperatures.  The  surface  of  sulfur-diffused  samples  generally 
showed  considerable  pitting  from  the  reaction  and  exhibited  appreciable 
losses  in  weight.  Again  data  were  obtained  by  the  incremental  sheet 
resistance  technique.  The  results  showed  that  the  diffusion  of  sulfur 
follows  a  complementary  error  function  distribution  in  gallium  arsenide. 

With  reservations  because  of  the  aforementioned  dlfficul ties>  the  data 
shown  plotted  in  Fig.  28  were  obtained.  One  value  on  this  curve  was 

1C 

obtained  by  the  radioactive  tracer  analysis  technique  of  S  .  This 
diffusion  coefficient  value  at  SSO"*  C  agrees  very  well  with  the  values 
obtained  from  the  sheet  resistance  technique. 
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Fig.  26 
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Fig*  27  Surface 


D(cm  /sec.) 


Fig.  28  Diffusion  of  sulfur  in  GaAs 
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H.  Magnesium  Diffusion  in  GaAs 


Diffusion  of  magnesium^  a  p«type  impurity^  in  gallium  arsenide 

was  studied  in  considerable  detail.  Diffusion  runs  were  made  at  800, 

900;  1000;  1100;  and  1150*’  C  in  sealed  ampoules.  The  amounts  of 

gallium  (placed  in  ampoule  to  prevent  copper  diffusion)  and  magnesium 

(40-60  atomic  percent;  respectively);  volume  of  ampoule;  and  size  of 

sample  were  kept  constant  throughout  the  runs.  Initially;  results 

were  evaluated  by  measuring  the  sheet  resistance;  p^;  of  the  p-layer 

and  the  Junction  depth;  x.  Values  of  the  surface  concentration;  C  ; 

were  obtained  from  a  set  of  curves  of  C  vs  p^x  for  various  starting 

o  s 

n-impurity  concentrations;  Cx.  These  curves  are  valid  if  the  impurity 

distribution  in  the  p-layer  is  assumed  to  be  a  complementary  error 

function  for  magnesium  diffusion  in  gallium  arsenide.  An  independent 

set  of  experiments  validated  this  assumption.  Figure  29  shows  a  plot 

of  C  vs  temperature  where  C  ranged  from  4.5  x  10^^  atoms/cc  at 
o  ift  ® 

800*’  C  to  5  X  10  atoms/cc  at  1150*’  C.  Calculation  of  the  diffusion 

-13  2 

coefficients;  D;  for  these  values  yielded  D  »  I  x  10  cm  /sec 
at  800*’  C  to  2.3  x  lO"^^  cm^/sec  at  1150*’  C.  A  plot  of  these  values 
of  D  vs  T  (temperature)  is  shown  in  Fig.  30.  Values  at  800*’  C  are 
believed  to  be  somewhat  in  error  because  of  inaccurate  measurement  of 
the  small  values  of  x  (<  I  micron  after  20  hours  of  diffusion).  There 
may  also  be  some  error  at  the  higher  temperatures,  where  considerable 
surface  loss  of  arsenic  Is  possible. 

Figures  31,  32  and  33  show  plots  of  concentration,  C.  vs  depth 
of  penetration,  x,  at  1100,  1000,  and  900“  C  for  magnesium  diffusion 
if  gallium  arsenide.  These  curves  tend  to  support  the  validity  of  the 
assumption  that  magnesium  obeys  a  complementary  error  function  distri¬ 
bution.  One  can  calculate  D  at  various  points  on  the  curves  and  find 
them  to  correlate  with  those  shown  in  Fig,  30.  The  results  In  Figs. 

31,  32  and  33  were  obtained  by  incremental  sheet  resistance  measurements, 
p^,  vs  depth  of  penetration,  x,  for  the  p-layers. 
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SURFACE  CONCENTRATION  (ATOMS/CC) 


1000/T(*K) 


TEMPERATURE  (‘C) 

Fig.  29  Diffusion  of  magnesium  in  GaAs  (n-type).  Source,  60  atomic  %  Mg. 
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DIFFUSION  COEFFICIENT  (CM^/SEC) 


1000/tCK)  - ^ 

M.P.  .7  .75  .80  .85  .90  .95 


TEMPERATURE  ('C) 

Fig.  30  Diffusion  of  magnesium  in  GaAs  (n-type).  Source  60  atomic  %Mg 
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CONCENTRATION  (ATOMS/CC) 


CONCENTRATION  (ATOMS/CC) 


DEPTH  IN  MICRONS 

Fig.  33  Concentration  at  a  function  of  diffusion  depth.  Source  iO  utomic  %Mg, 
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I.  Mercury  Diffusion  in  GaAs 


The  diffusion  of  mercury;  as  an  acceptor;  in  gallium  arsenide 

was  investigated.  Two  samples  of  gallium  arsenide  were  diffused  44 

203 

hours  at  1000'’  C  with  the  radioactive  Isotope  Hg  ;  at  a  mercury  pressure 
of  1,0  atmosphere.  Data  from  both  samples  are  shown  in  Fig.  34.  The 
limits  are  shown  for  a  90%  confidence  level  (see  Appendix  A);  within 

these  limits  the  diffusion  coefficient  may  be  no  greater  than  7*3  x 

*I4  2  203 

10  cm  /sec.  For  a  more  accurate  determination;  Hg  with  a  higher 

specific  activity  than  was  available  would  be  required.  The  value 
-14  2 

7.3  X  10  cm  /sec  was  determined  by  assuming  a  complementary  error 
function  distribution  and  finding  the  highest  0  that  would  fall  within 
the  predetermined  statistical  limits  of  the  data  taken.  This  is  not 
intended  as  an  actual  value;  but  merely  the  maximum.  Because  of  this 
low  diffusion  coefficient  for  mercury  at  1000^  C;  no  further  investiga¬ 
tion  of  mercury  diffusion  in  gallium  arsenide  was  deemed  necessary. 

J.  Thermal  Conversion 

A  study  of  contamination  as  a  cause  for  the  "thermal  conversion" 
of  gallium  arsenide  from  n-  to  p-type  during  diffusion;  was  accomplished. 
Copper  was  established  as  the  offending  agent.  The  exact  sources  of 
contamination  were  not  defined;  although  the  possibilities  are  numerous. 

In  the  experimental  procedure  a  slice  of  n-type  gallium 
arsenide  of  known  concentration  was  sealed  in  a  thoroughly  cleansed  quartz 
ampoule;  with  either  a  known  amount  of  gallium  or  arseniC;  at  a  vacuum 
of  about  one  micron.  In  each  run  the  preparation  of  slices  and  tubes 
was  identical;  samples  were  heated  at  1100^  C  for  one  hour  with  the 
same  cooling  cycle  on  each  run.  After  each  experiment  conductivity 
type  and  carrier  concentration  of  each  sample  were  ascertained  and  a 
spectrograph! c  analysis  made  for  existing  impurities.  In  the  case  of  a 
control  slice;  which  was  not  heated;  a  spectrographic  analysis  was  made  to 
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-IMPUNITY  CONCENTRATION  (ATOMS/CC) 


X-DEPTH  OF  PENETRATION  (MICRONS) 


Fig,  34  diffusion  in  GaAa 


determine  the  original  impurities  in  the  crystal.  Table  V  shows  the 
results  of  these  experiments^  from  which  the  following  conclusions 
can  be  made, 

(1)  Copper  (a)  is  not  found  in  the  grown  crystals  used  in 
the  experiments^  but  (b)  is  introduced  during  the  heat 
treatment  process. 

(2)  The  number  of  acceptors  added  during  the  treatment  is 
proportional  to  the  amount  of  copper  added. 

(3)  Various  chemical  treatments  of  the  gallium  arsenide 
surfaces  reduce  the  amount  of  copper  introduced^  but 
do  not  eliminate  it. 

(4)  Contamination  from  the  quartz  ampoule  appears  to  be  a 
larger  factor  in  sample  preparation, 

(5)  Excess  gallium  in  the  ampoule  significantly  reduces 
the  amount  of  copper  introduced  into  the  gallium  ar¬ 
senide  and  permits  the  use  of  n-type  material  with 
starting  concentration  as  low  as  5  x  10^^  atoms/cc  to 
make  transistors.  Placing  the  gallium  in  contact  with 
the  slices  does  not  aid  further  reduction. 

(6)  Excess  arsenic  in  the  system  has  much  less  effect 
than  gallium  on  copper  contamination. 

Attempts  to  eliminate  the  copper  contamination  completely 
have  not  yet  been  successful. 
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Control  slice  not  heated 
ND  stands  for 'hot  detectable" 

Cal  Hum  was  placed  in  actual  contact  with  the  sample 


IV.  SILICON  DIOXIDE  MASKING  ON  GALLIUM  ARSENIDE 


A,  General 

A  study  of  the  effect  of  reactively  sputtered  silicon  dioxide  films 
on  the  surface  of  gallium  arsenide  as  a  mask  against  impurity  diffusion  was 
initiated.  Preliminary  investigations  were  made  on  zinc  and  magnesium  dif- 

o 

fusion  through  12,000  A  layers  of  silicon  dioxide  on  partially  coated  n-type 
gallium  arsenide  slices.  Figure  35  shows  the  result  of  zinc  diffusion 
in  the  presence  of  the  film.  Zinc  was  diffused  from  a  1%  zinc-gallium 
alloy  in  a  sealed  ampoule  at  800'’  C  for  five  hours  and  showed  no 
penetration  on  the  protected  regions,  but  a  two-micron  penetration  on 
the  uncoated  regions.  Figure  36  shows  the  result  for  magnesium  dif¬ 
fusion.  The  magnesium  diffused  from  a  60%  magnesium-gallium  alloy 
in  a  sealed  ampoule  at  1000®  C  for  one  hour  showed  no  penetration  in 
the  area  of  the  silicon  dioxide  coating  but  penetrated  1.3  microns 
in  the  unprotected  region.  Diffusion  was  carried  out  for  these  impurities 

O 

with  thinner  films  (^s  2000  A),  through  which  visible  diffusion  was 
observed.  Hence,  it  can  be  concluded  that  the  masking  effect  is  de¬ 
pendent  on  the  silicon  dioxide  thickness  for  a  given  diffusion  tempera¬ 
ture  and  time.  To  fully  evaluate  this  silicon  dioxide  mask,  both 
n-  and  p-impurities  for  gallium  arsenide  should  be  thoroughly  investi¬ 
gated  by  the  radioactive  tracer  technique.  Results  should  indicate 
the  eventual  feasibility  of  a  gallium  arsenide  planar  transistor 
structure. 


From  this  initial  masking  study  it  was  observed  that  impurities 
may  be  effectively  diffused  through  this  mask  with  no  dissociation  of 
the  gallium  arsenide  during  the  course  of  the  diffusion.  With  gallium 
arsenide  and  other  compound  semiconductors  similar  in  volatility  at 
diffusion  temperatures,  the  surface  of  the  crystal  is  lost  in  varying 
degrees  whether  the  diffusion  is  accomplished  by  sealed  ampoule, 
semlclosed,  or  open-tube  technique.  To  depress  this  loss,  a  closed 
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Fig.  35  Zinc  diffusion  with  1%  Zn-Ga  alloy  in  sealed  ampoule 
at  800«C  for  5  hours .  Magnification  70X 
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Fig.  36  Mg  diffuaion  with  60%  Mg-Ga  alloy  in  sealed  ampoule 
at  1000«C  for  one  hour.  Magnification  140X 
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system  and/or  lower  diffusion  temperature  is  usually  used.  However, 
even  by  this  method  the  surface  loss  was  still  undesirable  for  electronic 
devices.  The  advent  of  this  oxide  mask  will  permit  two-zone  open- 
tube  diffusion  and  a  more  accurate  evaluation  of  the  semiclosed  system. 
Preliminary  investigations  in  both  cases,  in  the  presence  of  the  silicon 
dioxide,  showed  surface  loss  of  gallium  arsenide  to  be  nonexistent. 

In  a  simultaneous  diffusion,  two  slices  of  gallium  arsenide 
were  diffused  5  1/2  hours  at  800°  C  with  1%  Zn-Ga  alloy  in  the  semi¬ 
closed  system.  Both  slices  had  been  chemically  polished,  but  only 

0 

one  was  protected  with  a  2000  A  layer  of  silicon  dioxide.  After 
diffusion,  the  film  was  removed  with  hydrofluoric  acid.  The  surface 
under  the  oxide  layer  of  the  protected  slice  showed  no  evidence  of 
arsenic  loss.  The  unprotected  surface  was  covered  with  small  droplets 
of  gallium,  definitely  showing  loss  of  arsenic.  Junction  depths  measured 
on  both  slices  were  approximately  1.6  microns.  However,  diffusion 
depths  of  the  two  samples  cannot  be  directly  compared,  since  it  is 
impossible  to  determine  accurately  the  amount  of  surface  loss  from  the 
unprotected  slice. 

We  also  made  a  two-zone,  open- tube  diffusion  run  which 
further  verified  effectiveness  of  the  silicon  dioxide  as  protection 
against  surface  loss.  Figure  37  shows  the  results  very  dramatically. 

O 

The  slice  was  coated  in  stripes  with  a  12,000  A  layer.  Diffusion  was 
carried  out  wi'th  the  slice  at  800°  C,  zinc  source  at  360°  C,  and  a 
helium  carrier  gas  at  the  rate  of  1  cu  ft/hr.  No  visible  diffusion 
was  observed  in  either  the  protected  or  unprotected  areas.  However, 
the  surface  loss  in  the  uncoated  areas  was  12.5  microns,  with  none 
observable  in  the  coated  area.  These  results  open  a  whole  new  avenue 
of  diffusion  techniques  for  the  more  volatile  compound  semiconductors. 
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Fig.  37  Two- zona  open  tube  diffusion  at  800*^0 »  Zn  source  at 

360^C  in  He  flow  (ICHF)  for  1  hour  12,000®A  film  blocked 
diffusion 
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B,  Magnesium  Diffusion  Through  SIO^  Film  into  GaAs 

Magnesium  diffusion  into  gallium  arsenide  through  a  silicon 

dioxide  film  was  analyzed.  Initially^  an  experiment  was  performed  to 

evaluate  the  distribution  of  magnesium  in  both  the  film  and  the  gallium 

0 

arsenide.  The  slices  were  coated  with  a  1000  A  layer  of  silicon  dioxide 
on  chemically  polished  surfaces^  and  then  diffused  at  1000^  C  for  six 
hours  from  a  60%  Mg  ->  40%  Ga  source.  To  find  the  behavior  of  this 
diffusion^  the  concentration  in  the  film  was  determined  by  spectro- 
graphic  analysis  on  the  samples  prepared  by  dissolving  the  film  in¬ 
crementally  in  10%  HF.  In  the  gallium  arsenide  the  concentration  was. 
obtained  by  the  standard  Incremental  lapping  and  sheet  resistance 

- 12 

technique.  The  results  are  shown  in  Fig.  38,  where  a  0  •  6.76  x  10 
2 

cm  /sec  was  found  in  the  gallium  arsenide.  As  expected,  this  value  of 
D  is  in  fair  agreement  with  those  previously  reported  without  SiO^  films. 
This  film  should  have  no  effect  on  the  value  of  D  in  the  gallium  arsenide. 

Consequently,  a  series  of  diffusion  run'*  were  made  at  tempera¬ 
tures  from  800®  C  to  1150**  C,  and  diffusion  was  a^.ulyzed  on  both  pro¬ 
tected  and  unprotected  surfaces.  The  gallium  arsenide  slices  were  cut 
500  X  400  mils  from  approximately  0.01  ohm-cm  tin-doped  <lll>  grown 
crystals.  Mirror  finish  and  good  planarity  were  obtained  on  the  slices 
by  chemical  polishing,  and  the  possibility  of  impurity  contamination 
during  the  lapping  process  was  eliminated.  Thus,  final  sample  thickness 

o 

was  about  20  mils.  These  slices  were  then  coated  with  2000  A  thick 
silicon  dioxide  films  on  the  (T7T)  face  and  diffused  from  a  60%  Mg  - 
40%  Ga  source  in  an  ampoule  evacuated  to  approximately  one  micron.  The 
oxide  film  was  then  removed  with  HF,  and  effects  were  studied  by  measuring 
the  diffusion  parameters  on  both  surfaces.  Sheet  resistance  was  measured 
with  the  four-point  probe,  and  the  diffusion  depth  by  the  standard  lap- 
stain  technique.  The  surface  concentration  and  diffusion  coefficients 
were  calculated  assuming  the  complementary  error  function  distribution 
previously  shown  valid  for  magnesium  in  gallium  arsenide.  Figure  39 
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IMPURITY  CONCENTRATION  (ATOMS  PERC.C) 


DEPTH  IN  MICRONS 

Fig.  38  Oiffuaion  of  magneiium  in  GaAs  -  N-type  for  6  hours  at 
1000*C  (source  diluted  with  40%  gallium) 
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shoMS  the  results  of  this  investigation.  As  expected,  the  film  had  no 
effect  on  the  values  of  D.  Hovmver,  the  silicon  dioxide  definitely 
did  have  a  masking  effect;  Junction  depths  and  surface  concentration 
VMre  reduced  on  the  coated  surfaces. 

Appendix  B  of  this  report  presents  the  analytical  solution  for 
magnesium  diffusion  in  the  SlO^-GaAs  system. 

C.  Zinc  Diffusion  Through  SiO^  Film  into  GaAs 

Zinc  diffusion  through  a  SIO2  film  into  gallium  arsanlde  Mas 
studiad.  The  slices  of  gallium  arsenide  were  tin-doped  crystals 
to  a  concentration  of  1  x  10^^  atoms/cc  with  chemically  polished  surfaces. 

o 

Half  the  slices  were  then  coated  with  a  6500  A  thick  silicon  dioxide 
layer  and  one  coated  slice  and  an  uncoated  slice  placed  in  a  sealed 
ampoule  with  a  Zn-Ga  source.  In  all  cases  the  diffusion  temperature  of 
1000^  C  and  diffusion  time  of  18  hours  were  held  constant.  Only  the 
percentage  of  Zn  in  the  source  was  varied  in  the  three  runs  made.  The 
radioactive  Isotope  of  Zn^^  was  used>  allowing  tracer  technique  use 
for  the  diffusion  profile  analysis. 

Figures  40-42  show  the  concentration  profiles  for  a  1.0% 

Zn>  0.1%  Zn>  and  0.01%  Zn  source,  respectively.  In  all  cases  the  sur¬ 
face  concentrations  under  the  film  are  significantly  lower  than  on  the 
unprotected  surfaces.  The  depth  of  penetration  is  also  considerably 
less  on  the  protected  surfaces,  indicating  very  definitely  that  the 

SiO,  film  serves  as  a  mask  against  zinc  diffusion.  Mote  also  that 

^  18 
when  the  zinc  surface  concentration  is  less  than  10  atoms/cc,  the 

distribution  approximates  a  complementary  error  function.  Figures  4l 

and  42  show  diffusion  coefficients  calculated  for  these  cases,  and  they 

agree  very  well. 
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C-IMPURITY  CONCENTRATION  (ATOMS/CC) 


Fig.  40  Zinc  diffusion  in  GaAs 
-90- 


-IMPURITY  CONCENTRATION  (ATOMS /CC) 


-IMPURITY  CONCENTRATION  (ATOMS/CC) 


X-DEPTH  OF  PENETRATION  (MICRONS) 


Fig.  42  Zinc  diffusion  in  GaAs 
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The  validity  of  the  meager  data  exhibited  on  Fig.  42  is 
supported  by  electrical  measurements.  The  uncoated  slice  was  p-type 
on  the  surface  after  diffusion  but  the  masked  sample  was  still  n-type. 
ThuS/  in  regard  to  electrical  conductivity  this  sample  was  affectively 
masked  against  zinc;  but  because  the  tracer  technique  was  used,  zinc 
diffusion  at  a  level  below  the  bulk  donor  level  was  detected. 

Diffusion  coefficients  of  7  x  |0“*^  cm^/sec  and  6.5  x  lO"*^ 

2 

cm  /sec  shown  were  calculated  assuming  the  surface  concentration  under 
the  film  to  be  constant  during  the  time  of  diffusion.  Since  little  is 
known  about  the  actual  diffusion  of  zinc  through  the  Si02  film  itself, 
this  may  or  may  not  be  a  valid  assumption.  Work  will  be  carried  out 
to  determine  whether  the  surface  concentration  of  the  zinc  under  the 
film  is  time  dependent. 
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V.  MISCELLANEOUS  TECHNOLOGY  STUDIES  ON  GaAs 


Several  miscellaneous  technology  studies  made  on  gallium  arsenide 
included:  chemotaxial 1y  polishing  of  surfaces,  delineation  of  Junctions 
by  etch  or  stain  techniques,  post-alloy  diffusion  studies,  and  surface 
"passivation"  of  gallium  arsenide  transistors. 

A,  Chemotaxial  Polished  Surfaces 

Because  of  the  adverse  effects  of  surface  damage  to  device 
characteristics,  as  introduced  by  mechanical  polishing,  another  method 
for  preparing  gallium  arsenide  surfaces  was  sought.  In  particular, 
attention  was  focused  on  a  chemical  polishing  process  reported  by  the 
Bell  Telephone  Laboratories, 

The  process  which  finally  evolved,  consisted  of  placing  the 
gallium  arsenide  slices  between  two  teflon  discs  which  are  submerged 
in  a  solution  of  80  parts  H2S0^,  10  parts  H^O^  and  10  parts  H2O  in  a 
rotating  inclined  teflon  beaker.  Figure  43  is  a  picture  of  the  system 
as  presently  used. 

This  system  has  yielded  excellent  optically  polished  surfaces, 
free  of  surface  damage  and  negligible  lensing  on  the  outside  periphery 
of  the  slices.  An  important  observation  was  that  the  above  solution 
would  only  polish  the  (111)  or  arsenic  face  of  the  slice.  If  atten^ts 
were  made  to  polish  the  gallium  face,  a  very  rough,  pitted  surface  would 
result.  This  chemotaxial  polishing  technique  is  now  being  used  ex¬ 
clusively  to  prepare  slices  for  diffusion  and  device  studies. 

B.  Delineation  of  Junctions 

Many  etches  and  stains  were  investigated  for  delineating  the 
diffused  and  alloyed  junctions.  This  is  of  considerable  importance  In 
evaluating  diffusion  depths,  alloyed  regrowth  regions  and  the  base 
widths  in  transistors. 
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A  very  satisfactory  stain  for  decorating  diffused  p-n  Junctions 
In  gallium  arsenide  was  found.  It  consists  of  2  parts  HF,  2  parts 
acetic  acid,  4  parts  HNO^,  I  part  methanol,  and  I  part  H^O.  This 
mixture  is  more  effective  if  allowed  to  age  prior  to  use.  Wafers  con¬ 
taining  diffused  layers  are  lapped  and  polished  with  a  5"  bevel  and 
stained  to  reveal  the  junction.  The  nature  of  the  stain  is  such  that 
the  n-region  will  stain  black.  This  stain  is  now  used  in  evaluating 
all  diffusion  depths. 

An  etch  that  has  been  used  successfully  in  delineating 
Junctions  (both  diffused  and  alloyed)  consists  of  5  parts  HCI,  I  part 
HNOj,  and  2  parts  of  H2O.  This  solution  will  etch  the  p-region  pref¬ 
erentially.  Both  alloyed  p-n  Junctions  and  alloyed  emi tters-dif- 
fused  base  n-p-n  transistors  have  been  cross-sectioned  and  delineated 
with  this  etch.  Figure  44  shows  a  cross-section  through  the  active 
region  of  a  transistor  that  has  been  etched  with  the  above  solution. 

The  light  region  at  the  bottom  on  the  photograph  is  the  n-type  col¬ 
lector,  while  the  darker  band  is  the  diffused  p-type  base  region.  The 
light  stripe  penetrating  into  the  base  is  the  n-type  regrowth  region  which 
acts  as  the  emitter.  This  particular  transistor  had  an  emitter  diameter 
of  0.002  in.  with  a  base  width  of  approximately  3u.  This  process  is 
being  used  routinely  in  delineating  Junctions  in  diodes  and  transistors. 

C.  Post-Alloy  Diffusion 

Only  a  token  effort  was  applied  in  the  post-alloy  diffusion 
technique  in  gallium  arsenide.  Investigation  of  this  area  was  limited 
to  one  alloy,  Au-Sn-Zn.  A  major  problem  was  encountered  with  this 
alloy.  The  concentration  of  Zn  in  the  alloy  had  to  bo  maintained  at 
a  level  of  less  than  0.01%  to  avoid  compensation  of  the  Sn  in  the  n- 
regrowth  region.  This  is  because  of  the  2  orders  of  magnitude  difference 
In  the  segregation  coefficients  of  Zn  and  Sn.  The  net  result  is  that 
the  Zn  level  in  the  alloy  has  to  be  maintained  at  a  low  level  comparable 
to  other  unknown  and  probably  undesirable  impurities  present.  Conse- 
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Cross  section  of  the  active  region  of 
a  diffused  base-aiioyed  emitter  Galiium 
Arsenide  transistor.  Magnification  iOOO  X. 


quentlyj  very  erratic  results  are  obtained  with  this  alloy.  However^ 
several  transistors  which  exhibit  3's  ea  I  were  realized  by  this  technique. 
This  area  of  endeavor  was  replaced  by  more  promising  techniques. 


0.  Surface  "Passivation"  of  GaAs  Transistors 

A  srtwll  effort  was  made  to  improve  and  possibly  "passivate" 
the  state-of-the-art  gallium  arsenide  transistors  by  chemical  treat¬ 
ment.  As  is  known,  the  current  gain,  3,  of  the  transistors  falls  off 
rapidly  at  low  current  levels  and  is  very  surface  sensitive.  A  chemical 
etch  of  5  parts  NaOH  to  1  part  H2O2  improves  on  the  low  current  level 
3  and  Increases  the  value  by  as  much  as  an  order  of  magnitude.  Figure 
45  demonstrates  the  effect  of  a  20  second  etch  on  a  gallium  arsenide 
transistor.  The  low  current  level  3  is  seen  to  improve  dramatically 

after  this  etch,  with  3  increasing  from  100  to  480  at  5  volts  and  5  >na. 

This  improvement  may  be  ettributed  to  a  lowering  of  the  surface  recombina¬ 
tion  velocity  in  the  vicinity  of  the  emitter.  Device  stability  of  the 

units  was  also  improved  by  this  treatment.  It  was  observed  that  if 
more  than  a  20  second  etch  was  applied  no  further  Improvement  in  3  was 
realized,  but  that  the  saturation  voltage,  was  considerably  degraded. 
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Fig.  45 


1  ma/div 


VI.  TUNNEL  DIODES 


A.  General 


A  considerable  portion  of  the  device  effort  on  this  program 
was  devoted  to  the  evaluation  of  the  feasibility  of  GaAs  as  a  tunnel 
diode  material  and  development  thereof.  The  theory  of  operation  of 
tunnel  diodes  indicated  that  gallium  arsenide  should  be  superior  to 
both  germanium  and  silicon  in  all  respects.  Because  of  the  high  energy 
gap>  1.35  eV)  the  voltage  swing  over  the  negative  resistance  region 
should  be  greater  than  that  obtained  for  either  of  the  other  materials. 
Since  the  band  structure  of  gallium  arsenide  does  not  require  the 
introduction  of  phonon  interactions  to  conserve  momentum  during  tunneling 
transitions,  as  do  the  structures  of  both  silicon  and  germanium,  the 
probability  of  such  transistions  may  be  expected  to  be  higher  for  GaAs. 

In  addition,  the  probability  will  also  be  enhanced  by  the  lower  value 
of  the  electron  effective  mass  in  GaAs.  Therefore,  the  same  Junction 
dimensions  and  geometry  in  gallium  arsenide  should  result  in  higher 
peak  currents  and,  since  the  valley  current  does  not  depend  on  these 
factors,  should  yield  higher  peak-current  to  valley-current  ratios. 

As  is  known,  these  assumptions  were  verified  on  n-  on  p-type  gallium 
arsenide.  However,  as  Is  also  now  known,  a  serious  degradation  problem 
of  the  peak  current,  I^,  is  noted  on  these  units  at  elevated  temperatures 
and  operation  at  current  levels  of  greater  than  2  I^.  To  date  this 
problem  still  has  not  been  resolved  as  to  the  exact  mechanism  of 
failure  or  means  of  resolving  it. 

B.  Fabrication  Technology  and  Characteristics 

However,  since  gallium  arsenide  did  appear  to  be  a  favorable 
material,  attempts  were  made  to  construct  tunnel  diodes  from  It.  The 
experiments  involved  the  formation  of  sharp,  heavily  eloped,  p-n  junctions 
by  alloying  acceptor-containing  impurity  pellets  to  degenerate  n-type 
base  material  and  by  alloying  donor-containing  Impurity  pellets  to  de¬ 
generate  p-type  base  material. 
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Diodes  were  constructed  from  selenium  or  tel  lurium^-doped 
n-type  GaAs  by  alloying  impurity  elements  such  as  zinc  and  cadmium 
diluted  in  metals  such  as  gold,  silver  and  tin.  The  best  results  were 
obtained  using  the  selenium-*doped  ingot  and  alloying  a  90%  tin  -  10% 
zinc  pellet.  Peak- to- va I  ley  ratios  up  to  4:1  were  obtained.  The  fact 
that  the  tin-zinc  alloy  used  for  forming  the  heavily  doped  p-region 
yielded  the  best  results  is  somewhat  surprising  since  tin  itself  is  an 
excellent  donor  impurity.  It  is  probably  that  higher  peak-to-val ley 
ratios  were  not  obtained  with  p-  on  n-Junctions  because  the  doping  levels 
in  either  or  both  sides  were  not  sufficiently  degenerate. 

The  results  obtained  using  zinc-doped  p-type  base  material 

and  alloying  pellets  of  tin^  or  tin  containing  small  amounts  of  other 

donors  such  as  selenium  or  tellurium^  were  considerably  more  promising. 

The  best  characteristics  realized  were  on  p-type  gallium  arsenide  which 

19 

was  doped  to  a  level  of  about  5  x  10  zinc  atoms/cc.  A  degenerate  n- 
type  regrowth  region^  forming  a  sharp  Junction,  was  formed  by  alloying 
tin.  Generally  pellets  of  0.002  in.  nominal  diameter  were  alloyed  at 
a  temperature  of  550-600 C.  Rapid  cooling  was  used  to  limit  diffusion 
and  prevent  broadening  of  the  transition  region.  Diodes  made  by  this 
process  yielded  peak-to-val ley  current  ratios  as  high  as  70:1.  Figure 
46  shows  a  photograph  of  the  V-I  characteristics  of  one  such  diode* 

As  is  seen,  the  voltage  swings  in  the  negative  resistance  region  were 
generally  about  0.4  volts. 

One  problem  with  this  process  was  the  wide  variations  obtained  in 
peak  current  value.  This  was  and  still  is  a  inherent  problem  in  the 
fabrication  of  tunnel  diodes.  To  bring  these  peak  current  values  to  a  given 
level,  an  electrolytic  etching  technique  utilizing  a  5%  sodium  hydroxide  solu¬ 
tion  and  10  ma  etching  current  was  employed.  During  etching  the  peak  current 
value  was  monitored  and  when  the  desired  value  was  obtained,  etching  was  stopped. 
Table  VI  shows  typical  values  for  the  peak  current  after  fabrication  of  diodes 
and  after  etching  to  a  nominal  value  of  10  ma.  Generally,  peak  current 
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Fig.  46  V-I  characteristics  of  a  gallium  arsenide 
tunnel  diode 
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Table  VI 

Results  of  Controlled  Electrolytic  Etching 
of  Gallium  Arsenide  Tunnel  Diodes 


1 

P 

After  Etching 

t  No. 

(before  etchlnq) 

In  5%  NaOH 

1 

66 

9.7 

2 

290 

10.0 

3 

39 

12.0 

4 

50 

9.5 

5 

50 

II. 2 

6 

54 

9.6 

7 

46 

9.3 

8 

20 

8.0 

9 

8l 

12.0 

10 

25 

10.6 

II 

50 

9.0 

12 

50 

9.2 

13 

25 

10.2 

14 

29 

10.0 

15 

24 

10.0 
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densities  of  ebout  500  emperes/cm  were  obteined,  el  though  velues  of  en 
order  of  nugnltude  end  more  higher  had  been  realized.  Capacitancei 
measured  at  the  vaMey  point,  as  low  as  2-4  picofarads  for  a  peak  current 
of  10  me  were  obtained.  Switching  times  of  less  than  I  nanosecond  were 
also  realized. 

C.  Anomalous  Characteristics 

At  the  time  of  these  studies  several  anomalies  were  observed. 

On  some  diodes  a  hump  or  Inflection  wes  observed  on  the  V-I  curves  In 
the  diffusion  current  region.  Figure  47  shows  a  photograph  of  this 
effect.  Experiments  were  performed  In  an  attempt  to  find  the  source 
of  this  effect.  An  Ingot  of  degenerate  zinc-doped  gallium  arsenide  was 
purposely  saturated  with  oxygen  whan  prepared.  This  element  was  known 
to  cause  a  deep  lying  level  in  gallium  arsenide.  Tunnel  diodes  were 
constructed  by  the  normal  process.  On  a  few  of  the  units  a  hump  was 
again  observed  on  the  V-I  curve.  However,  here  It  occurred  near  the 
center  of  the  valley.  Figure  48  shows  a  photograph  of  such  a  curve. 

It  was  believed  that  these  humps  correspond  to  a  second  tunneling  process 
Into  e  deep  Impurity  or  defect  level  in  the  material.  Different  Im¬ 
purities  which  give  levels  in  different  positions  In  the  forbidden  gap 
would  causa  the  hump  to  occur  at  different  values  of  bias  voltage,  es 
was  observed,  A  degenerate  zinc-doped  sample  was  doped  with  copper  during 
the  compounding  operation  to  determine  whether  It  would  also  laad  to  a 
hump  In  the  characteristic  curve.  Although  no  such  affect  was  observed, 
diodes  constructed  from  the  material  had  peak  current  densities  from 
two  to  three  orders  of  magnitude  lower  than  normally  obtained.  This 
result  was  unexpectad  and  its  causes  were  not  understood. 

Another  Interesting  observation  concurrent  with  the  above 

phenomena  was  the  dacrease  In  peak/valley  current  ratio  with  Increasing 

zinc  concentration.  In  general,  the  degenerate  p-type  gallium  arsenide 

19 

used  for  the  diodes  had  a  zinc  concentration  of  about  5  x  10  /cc. 

If  this  concentration  was  Increased,  the  peak  current  on  fabricated 
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Fig.  47  V*I  Characteristics  of  a  Gallium  Arsenide  Tunnel  Diode  Shoeing 
Anomalous  Hump  in  the  Curve 
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Fig.  48 


V-I  Characteristics  of  a  GaHium  Arsenide  Tunnel  Diode  made 
from  Oxygen  Containing  Material  Showing  Anomalous  Hump  in 
the  Valley  Region 


diodes  increased  because  of  the  Increased  tunneling  probability.  Howevar, 
the  valley  current  also  increased  at  a  rate  sufficient  to  cause  a  decrease 
in  the  peak/valley  current  ratio.  This  increase  in  valley  current  was 
believed  to  result  from  an  effect  similar  to  that  caused  by  neutron 
irradiation,  where  the  same  behavior  was  reported.  Zinc,  in  these  very 
heavy  concentrations,  could  strain  the  crystal  lattice  sufficiently  to 
cause  the  introduction  of  defect  bends  in  the  forbidden  energy  gap. 
Additional  current  would  then  result  from  tunneling  into  these  defect 
bands,  the  carriers  subsequently  undergoing  a  transition  to  the  valence 
band.  In  addition,  some  interstitial  zinc  levels  could  be  introduced. 

This  explanation  for  the  aforementioned  phenomena  was  never  verified 
or  discredited. 

D.  Temperature  Behavior  of  Characteristics 

An  Interesting  study  that  we  made  was  to  evaluate  the  gallium 
arsenide  tunnel  diode  over  a  wide  temperature  range.  Although  the 
physical  properties  of  gallium  arsenide  should  allow  device  operation 
to  temperatures  in  excess  of  350*  C,  the  presence  of  the  previously 
mentioned  tin  placed  an  upper  limit  of  about  200*  C.  Hence,  considerable 
effort  was  spent  in  search  of  a  more  satisfactory  donor  alloy.  Sue- 
cessful  results  were  obtained  from  the  62%  Au  -  38%  Sn  (by  weight) 
alloy,  mentioned  under  alloy  studies  of  this  report.  However,  great 
care  had  to  be  taken  when  raising  these  devices  to  temperatures  in 
excess  of  300*  C.  Figure  49  shows  the  V*I  characteristics  for  a  GaAs 
tunnel  diode  for  temperatures  from  >196*  C  to  345*  C.  At  room  temperature 
the  unit  has  a  paak/velley  current  ratio  of  17:1  and  increases  to  32:1 
at  -196*  C.  At  34s*  C  the  ratio  decreased  to  about  4:1.  A  better 
understanding  of  this  can  be  obtained  by  looking  at  Fig.  SO  where  the 
peak  end  valley  currants  era  plotted  against  temperature.  The  valley 
current  is  seen  to  increase  monoatomical ly  with  an  increase  in  temperature, 
end  in  the  range  of  77*  K  to  300*  K  increases  by  SS%,  whereas  in  the 
same  temperature  range  the  peak  current  changes  by  only  10%.  Conse¬ 
quently,  the  peek- to- valley  ratio  increases  considerably  in  lowering 
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temperatures.  The  decrease  in  peak  current  at  lower  temperatures  Is 
expected^  since  the  energy  gap  Increases  and  hence  tunneling  decreases.' 
The  peak  current  decreases  at  higher  temperatures^  which  is  expected; 
because  of  the  Fermi  function  smearing  out  at  elevated  temperatures. 

The  degradation  problem  was  not  noted  in  this  unit;  probably  because 
this  unit  was  not  kept  at  elevated  temperatures  for  any  prolonged 
period. 
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VII.  TRANSISTORS 


A,  General 

The  ultimate  goal  of  this  program  was  to  prove  feasibility  of  and 
develop  a  gallium  arsenide  transistor.  In  particular^  feasibility  samples 
of  a  high  temperature^  low  power  transistor  and  a  high  temperature^  medium 
power  transistor  were  to  be  fabricated. 

This  interest  in  gallium  arsenide  transistors  was  generated 

by  some  of  the  superior  basic  electrical  properties  of  gallium  arsenide 

compared  with  germanium  or  silicon.  Some  important  transistor  parameters 

are  limited  by  these  basic  electrical  properties  of  the  specific  semi* 

conductor  used  for  a  device.  The  maximum  operating  temperature  depends 

primarily  on  the  concentration  of  thermally  excited  carriers^  a  function 

of  the  forbidden  energy  gap.  Germanium  and  silicon  transistors^  with 

energy  gaps  of  0.7  and  1.1  eV>  respectively^  are  operable  to  temperatures 

of  85^C  and  150”C>  respectively.  For  gallium  arsenide  the  energy  gap 

is  1.35  eV  or  a  maximum  device  operating  temperature  of  approximately 

400''C.  The  power  handling  capability  of  a  transistor  depends  on  the 

temperature  at  which  the  collector-base  Junctions  can  operate  and  is 

also  influenced  by  the  energy  gap.  Again^  based  on  electrical  properties^ 

gallium  arsenide  transistors  should  be  superior.  The  ability  of  a 

transistor  to  amplify  at  high  frequencies  depends  on  material  properties 

in  a  more  complex  manner;  however^  one  important  factor  is  the  time 

required  for  minority  carriers  injected  from  the  emitter  to  diffuse 

across  the  base  region.  This  carrier  diffusion  is  related  to  the 

mobility  of  carriers  in  the  semiconductor.  Germanium^  with  electron 

2 

((x  )  and  hole  (tx  )  mobilities  of  3900  and  1900  cm  /Volt-sac  has  yielded 
n  p 

transistors  operable  to  frequencies  in  excess  of  I  kMc.  Silicon,  with 
2  2 

u  -  1500  cm  /Volt-sec  and  u  -  500  cm  /Volt-sec,  is  limited  to  a  few 

n  p  2 

hundred  He.  Again  gallium  arsenide  with  p  rs  10,000  cm  /Volt-sec  and 
2  ^ 

ix  ««  450  cm  /Volt-sec,  should  be  superior  to  either  germanium  or  silicon. 
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B.  Fabrication  Techniques 


A  primary  difficulty  in  constructing  transistors  from  compound 
semiconductors  such  as  gallium  arsenide  results  from  the  need  of  extremely 
thin  base  regions  because  of  the  very  low  carrier  lifetimes.  Diffusion 
techniques  were  chosen  to  form  the  base  region  because  of  their  proven 
usefulness  for  forming  such  thin  regions.  Alloy  techniques  were  chosen 
to  form  the  emitter  region,  Also^  an  n*-p-n  structure  was  selected  to 
take  advantage  of  the  higiier  electron  mobility. 

Figure  51  shows  a  sketch  of  the  gallium  arsenide  transistor 

structure  that  yielded  feasibility  units.  The  starting  material  was 

Se^  Te^  or  Sn  doped  n*type  single  crystals  to  a  level  of  0.5  -  1  x  10^^ 

atoms/cm^  grown  by  the  Teal-Little  technique.  Wafers  from  these  crystals 

were  then  diffused  in  an  evacuated  ampoule  with  Hn,  Mg  or  Zn.  All  slices 

were  either  mechanically  or  chemically  polished  on  the  As  surface  prior 

to  diffusion.  This  yielded  Junction  depths  of  approximately  0.3  mils 

1 S  3 

and  a  surface  concentration  of  1  -  5  x  10  atom$/cm^.  The  p-type  layer 
on  the  unpolished  surface  was  completely  removed  by  lapping.  A  small 
mesa  was  formed  on  the  p-type  surface  by  masking,  using  the  Kodak  photo¬ 
resist  technique,  and  then  etching  in  CP-4.  These  mesas  were  either  of 
4  X  4  or  3  X  5  mi  I  geometry.  The  n-type  emitter  region  was  then  formed 
by  alloying  a  2  mi  1  Sn  pellet  to  within  2-3  microns  of  the  collector- 
base  junction.  Simultaneously,  a  99.8%  gold  -  0.2%  selenium  foil  on  a 
platinum  tab  was  alloyed  to  the  n-type  surface  of  the  wafer  to  form  the 
collector  contact.  This  operation  was  carried  out  at  about  550"C  in  a 
small  resistance  heated  furnace,  A  2  mi  I  diameter  wire  containing 
99%  gold  -  1%  zinc  was  then  alloyed  to  the  p-layer  to  form  the  ohmic 
base  contact.  The  units  were  then  given  a  brief  CP-4  etch  to  clean  the 
surfaces. 
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Fig. 51  Structural  view  of  gallium  arsenide 

npn  mesa  transistor 


Transistor  Characterization 


C. 


Using  the  above  process^  many  transistors  were  fabricated  and 
their  electrical  parameters  measured.  Figure  52  shows  the  common  base 
output  characteristics  of  a  gallium  arsenide  transistor.  The  collector 
current  is  I  ma/div  and  the  emitter  current  is  1  ma/step.  The  current 
gain>  ot,  is  O.985  at  a  collector  current  of  5  ma  and  collector  voltage 
of  5  Volts.  The  common  emitter  characteristics  of  the  same  transistor 
are  shown  in  Fig.  53*  Again>  the  collector  current  is  1  ma/div,  and  the 
base  current  is  20  ^a/step.  This  yields  a  current  gain,  of  60  at 
5  ma  collector  current  and  5  Volts  collector  Voltage.  Gains  as  high  as 
3000  have  been  observed  when  units  were  given  the  surface  "passivation" 
etch  described  earlier  in  this  report. 

Table  VII  shows  values  for  the  common  base  "h"  parameters,  at 
a  collector  base  voltage,  V*  »  5  Volts  and  emitter  current,  I  *  I  ma, 

C  D  C 

for  several  state*of-the*art  gallium  arsenide  transistors  with  a  small 
signal  of  1000  cycles.  The  current  gain,  h^^,  shown  in  the  first  column 
is  as  high  as  0.99.  Input  impedances,  h.^,  of  50  ohms  and  output 
admittances,  h^|^,  of  0.5  micromhos  were  generally  obtained.  The  last 
two  columns  show  values  of  matched  power  gain  calculated  from  the 
measured  "h"  parameters.  Common  emitter  gains  as  high  as  43.6  db  were 
obtained. 


Table  VIII  shows  coniplete  measurements  on  a  typical  gallium 

“8 

arsenide  transistor.  The  leakage  current,  was  2  x  10  amperes. 

L  OO 

This  value  of  is  somewhat  higher  than  generally  realized  for  a 

2  -9  -10 
IS  mil  collector  area,  where  currents  of  10  to  10  amperes  were 

more  commonly  observed.  The  collector-base  breakdown  voltage,  BV^0qi 

was  15  volts  and  collector-emitter  breakdown  voltages,  BV^^^ ,  15  volts. 

Collector  and  emitter  translation  capacitances  were  3*2  pf  and  1.8  pf, 

respectively.  The  saturation  voltage,  Vj.g,  is  seen  to  be  0,6  volts  et 

-  10  ma>  Ig-  1  ma.  Generally  the  better  state-of-the-art  units 


114- 


I 

-Vc=l  VOLT/div— ^ 

e  output  characteristics  of  a  GaAs 
^ «  ma/step 


-115- 


r»g.  S3  Common  emitter  output  characterietica  of  a  OaAa 
tranaistor  «  0.  02  ma/atep 
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Ic  =IOma  lB  =  lma  I.SVolts 

(Sat)  Ic  =IOma  lB  =  lma  O.GVolts 

VcE  =  5  Volts.lE  =5ma  10 


exhibited  V  "  0.5  to  1.0  Volt  at  the  above  mentioned  conditions.  This 
CE 

particular  unit  showed  an  hp^  value  of  10  at  *  5  ma  and  ■  5  Volts. 

D,  *^h**  Parameters  vs  Frequency  Characterization 

To  adjust  transistor  fabrication  it  is  necessary  to  know  which 
parameters  are  limiting  transistor  performance.  This  type  information 
can  be  obtained  readily  from  electrical  measurements  of  the  transistor 
impedance^  voltage  transfer  ratio>  and  current  transfer  ratio  which  are 
related  to  the  physical  parameters  of  the  device.  For  example^  from  a 
plot  of  hj^  vs  frequency  an  accurate  value  of  base  resistance,  r|^  is 
obtained. 


Figures  5^  and  55  show  plots  of  the  real  and  imaginary  parts 
of  h.  vs  frequency  for  two  state-of-the-art  gallium  arsenide  transistors. 
The  result  of  such  plots  should  yield  a  simicircular  graph,  as  seen  in 
the  figures.  In  these  plots  the  lower  real  resistance,  the  point  at 
which  the  imaginary  component  becomes  zero,  is  a  good  measure  of  r^  *1*  r^. 

Since  r  is  the  emitter  series  resistance  and  is  usually  less  than  a  few 

e 

ohms,  the  measured  value  closely  approximates  r^.  Consequently,  values 
of  r^  shown  in  the  graphs  are  approximately  20  and  30  ohms,  respectively. 
These  were  not  unusual  units;  many  other  units  that  were  measured  ex¬ 
hibited  values  in  the  range  of  20  to  60  ohms.  These  measurements  were 
made  on  a  Transfer  Function  Bridge  for  the  high  frequency  end  and  on  an 
RX  meter  in  the  low  frequency  range. 

These  values  differ  by  almost  an  order  of  magnitude  from  the 
numbers  obtained  from  r^  C  measurements,  In  which  rf  fin  100  to  1000  ohms, 

DC  D 

The  measurements  obtained  from  h.  vs  frequency  are  believed  to  be  more 

1  e 

realistic.  Assuming  that  the  component  of  r'  under  the  emitter  is  low, 

D 

then  certainly  the  surface  resistance  between  the  emitter  and  base  should 

18 

b«  low  becauf*  of  the  high  surface  concentrations^  of  I  to  5  x  10  /cc. 
The  discrepancy  between  the  two  r^  values  is  not  understood  and  hence  has 
not  been  resolved  yet. 
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Fig.  54  Frequency  dependence  of  grounded  emitter 
input  impedance  for  on  experimental  gallium 
arsenide  transistor  V(sg  =  3V,  5ma. 
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Fig.  55  Frequency  dependence  of  grounded  emitter 
input  impedance  for  an  experimental  gallium 
arsenide  transistor  Vro  =  3V,  I-  5ma  . 


Figure  56  shows  a  plot  of  common  emitter  current  gain/ 
vs  frequency  for  a  state-of-the-art  gallium  arsenide  transistor.  If 
such  a  curve  falls  off  at  the  rate  of  6  db  per  octave  over  a  given 
frequency  range/  it  yields  a  value  for  the  gain-bandwldth-product/  f^. 
f^  is  the  frequency  at  which  the  value  of  h^^  *  0  db/  or  3  «  I,  As 
seen  from  Fig.  56/  the  low  frequency  value/  20  db/  corresponds  to  a  3 
of  10.  Beyond  100  Me  the  gain  does  fall  off  with  the  expected  6  db/octave 
slope/  thus  resulting  in  an  f^  value  of  approximately  730  Me.  In  turn  it 
can  be  safely  stated  that  the  or-cutoff  frequency  for  this  transistor  is 
for  ^  I  kMc.  f^  is  the  frequency  at  which  the  common  base  current  gain  is 
0.707  of  its  low  frequency  value  when  o  ^  I.  A  number  of  units  have  been 
measured  that  exhibited  f^  values  ranging  from  600  to  8OO  Me  as  shown 
in  Table  IX.  These  are  the  highest  values  of  f^  known  to  have  been  ob¬ 
served  on  gallium  arsenide  transistors  to  this  date. 

The  f^  of  a  graded  base  transistor  is  related  to  the  carrier 
diffusion  constant  and  the  transistor  base  width  by  the  approximate 
expression  ^  ^ 

(32) 

where  the  factor  of  2  is  an  enhancement  factor  due  to  the  graded  base.  Sub¬ 
stituting  values  obtained  on  experimental  transistors,  f^  ■  1000  Me  and  W  ■ 
3)1,  which  were  obtained  by  delineation  techniques,  gives  0^  115  cm^/sec. 

This  yields  an  average  electron  mobility  in  the  p-type  base  region  of 

u  ■  0  ■  4600  cm^/Volt-sec,  (33) 

n  Ki  n 

which  is  approximately  the  value  obtained  from  Hall  measuramonts  on  some 
crystals  used  in  fabricating  the  higher  frequency  units.  These  observed 
electrical  and  physical  parameters,  along  with  common  base  current  gain, 
or,  allow  an  estimate  of  the  electron  lifetime  in  the  base  region.  If  we 
assume  that  the  currant  gain  is  limited  by  base  transport.  It  is  given 
by 
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where  L  is  the  diffusion  length  of  electrons  in  the  base  region.  For 
a  ■  0.9  and  W  >  L  ■  6.7  x  10  cm  is  obtained.  The  diffusion  length 
may  be  expressed  as 

L  -  /D?  .  (35) 

Substituting  the  calculated  values  for  L  and  D  yields  a  lifetime 

-9  ^ 

T  ^  5  X  10  sec.  This  value  of  t  is  about  an  order  of  magnitude  lo^er 

than  what  has  been  calculated  by  the  same  technique  on  some  earlier  units 

with  an  f^  ^  300  me.  Nevertheless^  it  is  consistent  with  what  is  generally 

believed  to  be  the  range  of  carrier  lifetime  in  III-V  compound  materials. 

However^  this  difference  in  values  of  t  may  not  be  truly  realistic  in  that 

emitter  efficiencies  of  the  two  units  could  be  vastly  different.  It  has 

not  been  determined  what  the  emitter  efficiencies  are  in  gallium  arsenide 

transistors  and  the  assumption  that  y  1  may  be  erronous^  where 

Of  *  cr*  P  Y  >  (36) 

Pc  Pn 

ork  m  ^  — — C-)  in  the  collector  efficiency^ 

”c  ^n 

w2 

P  “  (1 - 5“)  the  base  transport  factor, 

■-nb 

P  W 

and  Y  *  (^  +  ^  i — )  the  emitter  injection  efficiency. 

^b  pe 

If  Y  not  really  approach  unity  very  closely,  as  assumed,  then  the 
original  assumption  that 

or  -  P 

is  in  error  and  the  value  of  t  calculated  is  meaningless.  Questions  have 
been  raised  regarding  emitter  injection  efficiency  on  alloyed  Junctions 
in  gallium  arsenide  transistors.  This  could  explain  the  discrepancy  in 
the  values  of  t  calculated  from  two  different  transistors. 
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It  would  be  very  desirable  to  be  able  to  determine  t  in  gallium 

arsenide  by  some  other  technique.  To  date  carrier  lifetime  has  not  been 

measured  meaningfully  by  any  of  the  commonly  known  techniques.  Photo- 

conductive  and  photomagnetoelectric  methods  have  been  attempted^  but 

neither  can  measure  the  apparent  short  lifetime  of  gallium  arsenide 
—8 

10  sec).  The  PME  approach^  however^  appears  more  promising.  Earlier 

difficulties  with  this  method  were  caused  by  the  inhomogeneity  of  gallium 

arsenide  and  the  extremely  high  surface  recombination  velocity  thought  to 

be  present.  Recently^  both  problems  are  believed  to  have  been  reduced 

and  the  PHE  method  will  be  tried  again.  As  pointed  out^  transistor 

current  gain  on  some  gallium  arsenide  units  would  seem  to  indicate 
-8 

T  10  sec;  however,  carrier  lifetime  is  probably  much  less  in  the 
majority  of  the  crystals,  since  good  3  units  are  generally  not  realized. 
Oxygen  content  of  the  material  is  known  to  introduce  traps  in  the  center 
of  the  forbidden  energy  gaps  of  gallium  arsenide  and  could  be  a  lifetime 
killer.  This  impurity  could  vary  greatly  from  crystal  to  crystal. 

Attempts  are  now  being  made  to  resolve  this  lifetime  problem  and/or 
emitter  injection  efficiency  in  gallium  arsenide. 

2 

The  previously  calculated  value  of  carrier  mobility,  \i  »  4600  cm  / 

volt-sec,  in  the  diffused  p-layer  is  of  interest,  since  it  points  up  the 

potential  frequency  capabilities  of  gallium  arsenide  as  a  transistor 

material.  This  value  is  a  factor  of  three  greater  than  calculated  for 

the  earlier  unit  that  exhibited  an  f^  ■  300  Me.  This  is  not  surprising 

in  that  the  later  unit  also  is  a  factor  of  three  better  in  f  for 

.  cr 

approximately  the  same  base  width  value,  W,  of  3  microns.  Consequently, 

2 

if  it  were  possible  to  attain  the  theoretical  limit  of  at  10,000  cm  / 
volt-sec  in  the  p-layers  and  W  at  0.25  -  0.50  microns,  as  in  the 
presently  highest  frequency  germanium  units,  gallium  arsenide  transistors 
would  far  surpass  any  presently  available  transistors  in  frequency. 

The  problem  remains  in  technology,  which  on  gallium  arsenide  has  not 
progressed  far  enough  to  allow  the  fullest  exploitation  of  GaAs  capability. 
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E. 


Switching  Time  Characterization 


Because  of  the  extremely  low  minority  carrier  Ufetimej  j,  ex¬ 
hibited  by  gallium  arsenide^  generally  two  orders  of  magnitude  lower  than 

<  »8  -6 

in  germanium  or  silicon  (  10  sec  compared  to  ^  10  sec)^  considerable 

interest  existed  in  the  switching  capabilities  of  the  state-of-the-art 
transistors.  Primarily^  considerable  reduction  in  the  storage  time^  t  ^ 
which  is  dependent  on  t  in  the  collector  region,  should  have  been  realized 
in  the  gallium  arsenide  transistor  vs  non-epi taxial  germanium  or  silicon 
transistors.  This  reduction  was  realized,  as  wi 1 1  be  shown  in  the  data 
to  follow;  however,  the  advantage  was  short-lived  with  the  advent  of 
epitaxial  germanium  and  silicon  transistors. 

Table  X  compares  the  switching  characteristics  of  some  state- 
of-the-art  units  with  those  of  five  2N706  silicon  mesa  transistors,  which 
were  designed  specifically  for  high  speed  switching  applications  at  that 
time.  These  values  were  measured  in  the  saturated  using  the  circuit 

shown  in  Fig.  57-  The  gallium  arsenide  units  were  significantly  faster 
than  the  2N706's.  Most  gain  was  in  the  storage  time,  where  the  gallium 
arsenide  units  had  values  between  1  and  10  nanoseconds  compared  with  28 
to  3^  nanoseconds  for  the  silicon  transistors.  The  delay  time,  t^,  and 
fall  time,  t^,  for  gallium  arsenide  were  also  somewhat  shorter;  while 
the  rise  time,  t^,  was  comparable  to  the  other  two  materials.  When  the 
switching  conditions  in  the  circuit  were  optimized,  a  further  reduction 
in  total  switching  time,  was  realized  in  Table  XI. 

As  the  geometry  of  the  state-of-the-art  gallium  arsenide 

transistors  was  reduced,  a  further  reduction  in  t^  ^  ,  was  realized. 

tota  I 

Several  of  these  units  were  measured  in  the  saturated  rTX>de  on  a  new 
switching  time  set  assembled  under  Contract  NOrd  18902  and  shown  in 
Fig.  58.  This  circuit  was  built  on  a  silver-plated  chassis,  with  highest 
quality  components  and  the  best  high  frequency  wiring  techniques;  pure  sil¬ 
ver  tab  leads  were  used  throughout.  This  circuit  was  used  with  the  sampling 
oscilloscope  and  a  fast  rise  time  pulse  generator,  shown  in  the  block 
diagram  of  Fig.  59.  With  considerable  effort  to  achieve  proper  resistive 
terminations  and  to  eliminate  strong  capacitances,  the  circuit  response 
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32 


Fig.  59  Block  diagram  of  the  test  equipment  used  to  observe  GaAs  transistor 
sw  itching  characteristics 


was  brought  within  the  rise  time  capability  of  the  oscilloscope  (0.8  nsec 

or  8  X  10~^^  sec).  One  gallium  arsenide  unit;  measured  in  this  test; 

set;  showed  a  total  switching  time  of  10.3  nanoseconds.  Figure  SO 

shows  the  collector  current  turn  on  response  for  this  unit  to  be  t^^  • 

6.5  nsec;  Fig.  61  shows  the  turn-off  response  as  t^^^  -  3.8  nsec,  giving 

the  t  ,  quoted  above.  The  current;  I..;  switched  was  10  ma  with  a  turn- 

on  current,  Ig|,  of  6.1  ma  and  a  turn-off  current,  1^2^  of  9.3  ma.  The 

turn-on  and  off  currents  were  optimized  to  produce  the  best  rise  and 

fall  times  for  the  unit  -  t^  =  5.2  nsec  and  t^  »  2.4  nsec.  This  value 

of  t  ,  was  the  best  realized  on  any  gallium  arsenide  units, 
total  ® 

As  was  expected  the  values  of  storage  time,  t^,  were  very  low, 

on  the  order  of  I  nanosecond  or  less.  Further  improvement  of  t^  and  t^ 

could  have  probably  been  realized  with  further  refinement  of  geometry. 

However,  the  t^  ^  .  would  still  have  been  of  about  the  same  order  of 
'  total 

^total'  2-5  nsecs,  which  is  now  being  realized  on  good  germanium  epitaxial 
transistors.  Consequently,  the  switching  advantage  that  was  believed 
to  be  gained  from  gallium  arsenide  transistors  vanished  and  work  on  this 
facet  of  the  transistor  program  was  discontinued. 

F.  High  Temperature  Characterization 

One  of  the  most  important  advantages  to  be  gained  from  gallium 
arsenide  as  a  transistor  material  is  its  high  temperature  capability  to 
operate  at  temperatures  up  to  400*'C.  All  these  units  were  built  with 
tin  alloyed  emitters,  which  limited  the  upper  limit  of  operation  to 
about  200"C,  Just  slightly  below  the  melting  point  of  tin.  Consequently, 
an  alloy-bearing  tin  was  sought  that  would  allow  higher  temperature 
operation.  It  was  believed  that  the  62%  Au-  38%  Sn  by  weight  alloy 
(melting  point  4l8"C)  would  resolve  this  problem.  However,  it  was 
found  that  upon  alloying  this  Au-Sn  compound  with  gallium  arsenide  the 
remaining  Au-Sn  became  liquid  at  about  SOO'^C,  Only  with  considerable 
care  in  handling  of  units  made  with  this  alloy  was  it  possible  to  obtain 
some  data  on  transistors  at  temperatures  up  to  350**C. 

Beta  vs  temperature  data  up  to  350'^C  was  obtained  on  a  number 
of  state-of-the-art  gallium  arsenide  units  made  with  the  Au-Sn  alloyed 
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t  s  2  nanoseconds /div 


Fig.  60  A  gallium  arsenide  transistor  collector 
current  turn  on  response,  to^*  6.5  nsec. 
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t  s  2  nanoseconds /div 


Fig.  61  A  gallium  arsenide  transistor  collector 
current  turnoff  response,  toff  ■  3.8  nsec. 
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emitter.  An  initial  determination  on  a  canned  transistor  at  temperatures 
up  to  350*C  showed  considerable  beta  degradation  after  returning  to  room 
temperature.  It  was  impossible  to  determine  whether  the  source  of  de¬ 
gradation  was  a  bulk  or  surface  phenomenon,  since  the  unit  could  not  be 
opened  for  study.  Hence,  further  investigations  were  carried  out  in 
vacuum.  Again  the  degradation  was  noted;  however,  each  time  the  beta 
could  be  returned  to  its  original  value  by  etching  (5  parts  NaOH-1  part 
H2O2).  This  condusively  showed  that  the  degradation  phenomenon  was  not 
in  the  bulk  material  and  was  not  permanent,  as  observed  in  the  gallium 
arsenide  tunnel  diode.  Figure  62  shows  the  temperature  effect  on  the 
common  emitter  output  characteristics  for  one  unit  while  cycling  from 
25*C  to  350*C  and  back  to  35*0.  As  noted,  very  little  change  took 
place.  Figure  63  shows  an  actual  plot  of  the  dc  beta  vs  temperature  for 

the  same  transistor  at  temperatures  up  to  SSO’C  (v  ■  5  volts  and  I  ■ 

C6  c 

2  ma).  As  noted^  beta  changed  very  little  over  this  temperature  range^ 
with  a  change  from  a  ^  ^  \k  to  ^  11.  This  probably  indicates  that  the 
minority  carrier  lifetime  over  this  temperature  range  is  relatively  un¬ 
changed.  Also  an  interesting  conclusion  is  that  the  frequency  over  this 
range  will  be  relatively  unchanged.  Independer.t  measurements  of  over 
this  temperature  range  in  bulk  material  showed  decreases  by  a  factor  of 

two.  Hence^  since  in  Eq.  (32)  f  «  0  ,  and  in  Eq.  (33)  D  «  Tp  it  is 

on  n  n 

apparent  that  0^  would  be  unchanged  from  25°C  to  350*0;  therefore,  the 
value  of  f^  should  not  change.  An  interesting  experiment  to  substantiate 
this  would  be  a  frequency  vs  temperature  study.  This  and  other  parameters 
vs  temperature  are  planned  for  gallium  arsenide  transistors. 

G.  Rellabi II  tv 

Twenty-one  state-of-the-art  gallium  arsenide  units  were  placed 
on  life  evaluation  tests  to  determine  the  kind  of  failure  that  may  be 
expected.  They  were  subjected  to  a  storage  test  at  I50*C  and  were  checked 
at  250  and  1000  hrs.  The  parameters  evaluated  were  at  5.0 

volts;  beta  at  1^  -  5  ma,  ■  5  Volts;  BV^^q  at  -  100  pa  and 
BVgfio  "  100  pa.  Storage  was  limited  to  I50*C  because  the  tin 

alloyed  emitters  on  the  units  would  melt  at  232*C.  These  units  were 
arbitrarily  chosen  for  the  test;  no  special  precautions,  such  as  canning 
in  drying  agents  or  inert  atmospheres  were  taken  in  their  preparation. 
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Tables  XII>  XIII^  XIV^  and  XV  show  the  effect  of  this  tempera¬ 
ture  storage  on  iQ^Qf  beta,  BV^0q>  ^nd  respectively.  Three  units 

failed  catastrophically  in  respect  to  3  and  after  1000  hrs,  with 

emitter-base  shorts.  Table  12  shows  that  the  values  in  general 
decreased  after  1000  hrs  storage.  This  might  be  from  drying.  However, 
in  many  cases  the  units  showed  no  further  decrease  beyond  250  hours,  but 
rather  a  slight  increase.  Perhaps  a  drying  agnnt  in  the  cans  would 
prevent  this.  Only  three  units  showed  a  very  large  change  in  beta  with 
the  bake-out.  Each  time  the  beta  decreased,  probably  because  of  some 
contaminant  on  the  device.  BV^^q  and  BV^^^  exhibited  very  little  change 
during  the  drying  process. 

The  primary  conclusion  to  be  drawn  from  these  results  in  that 
I50*^C  storage  has  no  drastic  effect  on  gallium  arsenide  transistors. 
Again,  certainly  no  real  degradation  was  observed  at  this  temperature, 
in  contrast  to  the  case  of  the  gallium  arsenide  tunnel  diode.  Not  much 
degradation  was  observed  in  beta  at  this  temperature,  as  at  350'’C.  How¬ 
ever,  in  either  case  the  changes  noted  were  conclusively  shown  to  be 
surface  phenomena. 
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TABLE  XII 


150*C  LIIE  STORAGE  EFFECT  ON  I, 


IcBo  Microamps  at  ■  5.0  Volts 

Unit 

0  Hours 

250  Hours 

1,000  Hours 

1 

0.028 

0.023 

0.019 

2 

0.006 

0.019 

0.04l 

3 

0.011 

0.00042 

0.0048 

4 

0.020 

0.010 

0.021 

5 

0.030 

0.0006 

0.012 

6 

0.045^ 

0.0015 

0.612 

7 

0.004 

0.00073 

0.0034 

8 

5.20 

4.21 

4.29 

9 

3.70 

1.05 

1.00 

10 

0.240 

0.200 

0.239 

11 

0.002 

0.004 

0.0075 

12 

152.00 

115.00 

129.00 

13 

36.30 

154.00 

153.00 

14 

0.001 

0.00067 

0.00079 

15 

0.01 

0.01 

0.0093 

16 

0.00076 

0.00063 

0.00034 

17 

0,00083 

0.00043 

0.00038 

18 

0.0036 

0.0020 

0.00145 

19 

0.0006 

0.00023 

0.00043 

20 

0.040 

0.0043 

0.0061 

21 

0.030 

0.020 

0.025 

)40 


TABLE  XIII 


150*C  LIFE  STORAGE  EFFECT  ON  BETA 


Beta  at  ■ 

5  ma,  Vj,j,  ■  5.0  volts 

Unit 

0  Hours 

250  Hours 

IjOOO  Hours 

1 

1.3 

1.4 

1.3 

2 

4.1 

2.6 

2.6 

3 

1.8 

1.8 

1.7 

k 

7.1 

6.2 

6.2 

5 

2.0 

2.0 

2.0 

6 

38.0 

4l.O 

45.0 

7 

3.8 

3.6 

3.6 

8 

7.1 

6.2 

6.2 

9 

1.8 

1.4 

1.4 

10 

2.0 

1.8 

2.0 

11 

7.7 

6.7 

7.1 

12 

1.3 

1.3 

0 

13 

4.6 

.9 

.7 

l4 

3.8 

3.3 

3.1 

15 

71.0 

71.0 

71.0 

16 

28.0 

31.0 

31.0 

17 

2.0 

4.1 

0 

18 

4.5 

1.8 

1.9 

19 

50 .0 

50.0 

50.0 

20 

17.0 

13.0 

0 

21 

46.0 

18.0 

18.0 
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TABLE  XIV 


150*  C  LIPE  STORAGE  EFFECT  ON  BV, 


BVggQ  in  Volts  at  -  100  jia 


Unit 

0  Hours 

250  Hours 

1.000  Hours 

1 

14.0 

14.5 

14.4 

2 

11.2 

11.2 

11.2 

3 

14.2 

14.6 

14.5 

4 

13.2 

13.6 

13.4 

5 

13.8 

14.0 

13.7 

6 

10.8 

11.2 

11.0 

7 

10.6 

10.5 

10.5 

8 

15.2 

15.9 

15.9 

9 

12.0 

11.0 

11.0 

10 

13.2 

13.1 

13.1 

11 

12.0 

12.5 

12.4 

12 

3.6 

4.0 

3.7 

13 

6.2 

4.4 

3.6 

14 

12.4 

12.8 

12.7 

15 

11.9 

12.0 

11.8 

16 

12.5 

13.2 

13.1 

17 

.  14.0 

14.6 

14.5 

18 

12.0 

12.5 

12.3 

19 

11.8 

12.2 

12.0 

20 

12.3 

12.5 

12.4 

21 

22.2 

22.1 

22.2 
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TABLE  XV 

150*C  LIFE  STORAGE  EFFECT  ON  BV_ 


BVgBo  in  Volts  at  -  100  ua 

Unit 

0  Hours 

250  Hours 

1,000  Hours 

1 

3.3 

4.0 

4.0 

2 

3.0 

1.1 

2.0 

3 

4.1 

3.9 

3.9 

4 

6.4 

6.4 

6.4 

5 

5.3 

4.5 

5.4 

6 

3.4 

3.2 

3.4 

7 

2.1 

3.0 

3.1 

8 

5.0 

4.5 

4.9 

9 

4.5 

3.6 

3.3 

10 

3.8 

4.8 

5.2 

11 

4.4 

4.1 

4.5 

12 

1.9 

2.5 

0 

13 

4.3 

2.3 

2.2 

l4 

3.9 

3.5 

3.7 

15 

6.5 

6.6 

6.7 

16 

5.5 

5.8 

6.0 

17 

6.1 

6.2 

0 

18 

6.3 

6.4 

■  6.4 

19 

-  ■ 

6.3 

6.4 

20 

5.0 

4.4 

5.7 

21 

4.9 

5.1 

5.2 
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VI It. CONCLUSION 


Texas  Instruments  successfully  executed  Part  I  of  the  program  as 
outlined  in  the  Foreword.  Alloy  contact  studies  for  both  ohmic  and 
rectifying  contacts,  diffusion  studies  of  primarily  acceptors,  and 
etching  studies  were  carried  out.  They  permitted  the  fabrication  of 
feasibility  transistors.  Some  typical  parameters  realized  on  these 
transistors  were  those  show  in  Table  VIII. 

The  requirements  of  Parts  2  and  3  of  the  program  were  only  partially 
fulfilled  because  no  satisfactory  high  temperature  and  high  injection 
efficiency  emitter  alloy  has  been  found.  Primarily  both  tin  and  gold-tin 
make  good  high  injection  efficiency  emitter  Junctions,  but  neither  allow 
operation  of  the  devices  up  to  the  400 **C.  However,  other  program  re¬ 
quirements  such  as  frequency,  switching  times,  and  power  gain  were  met 
and  measurements  of  the  feasibility  units  were  made,  h^^  vs  frequency 
measurements  at  V^g  of  3  Volts  and  1^  of  5  ma  yielded  gain  bandwidths 
products,  fj,  as  great  as  800  He.  This  value  is  the  highest  known  to 
have  been  observed  on  gallium  arsenide  transistors  and  is  believed 
limited  only  by  geometry,  h  parameter  measurements  at  f  of  1000  cps, 

Vqq  of  5  volts  and  1^  of  I  ma  yielded  common  emitter  power  gains,  6, 
ranging  from  30-43.6  db.  The  switching  characteristics  of  a  number  of 
these  units  were  also  investigated.  In  the  saturated  mode,  with  no  speed¬ 
up  capacitor,  a  number  of  GaAs  transistors  yielded  total  switching  times, 
ttotai'  l^~30  nanoseconds.  Host  of  this  time  was  in  the  rise  time, 

t^,  and  fall  time,  t^,  and  not  in  storage  time,  t^,  which  was  typically 
less  than  one.  Several  units  with  gold-tin  alloyed  emitters  were  heated 
to  350’C  with  very  careful  handling.  The  betas  of  these  transistors 
were  found  to  change  negligibly  over  the  temperature  range  of  25-350*0. 
However,  frequency  vs  temperature  over  this  same  range  could  be  expected 
to  decrease  by  a  factor  of  two. 

In  essence,  this  program  has  demonstrated  the  feasibility  of  gallium 
arsenide  transistors.  It  has  shown  the  tremendous  potential  of  GaAs 
devices  to  surpass  either  silicon  or  germanium  devices  for  higher  tempera¬ 
ture  operation,  higher  power  capabilities,  and  increased  frequency 
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response.  However,  before  this  potentiel  can  be  realized,  edvances  in 
GaAs  materials  and  device  technology  comparable  with  those  in  silicon 
end  germanium  device  technology  wi 1 1  be  necessary. 
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APPENDIX  A 

ESTIMATING  STATISTICAL  ERRORS  IN  DETERMINING 
NUCLEAR  DISINTEGRATION  RATE 


ESTIMATING  STATISTICAL  ERRORS  IN  DETERMINING 


NUCLEAR  DISINTEGRATION  RATE 


In  studying  diffusion  in  solids  using  radioactive  elements,  a  measure  of 
the  rate  of  nuclear  disintegrations  in  samples  of  known  sizes  is  obtained  to 
determine  the  concentration  of  the  tracer  element*  Since  the  decay  of  radio¬ 
active  elements  takes  place  in  a  random  manner,  errors  are  introduced  in  count¬ 
ing  rate  determinations  by  the  random  variations.  Consideration  of  the 
uncertainty  in  the  counting  rate  from  random  fluctuations  is  particularly 
important  for  low  activity  samples,  because  the  likely  statistical  errors  may 
far  exceed  the  errors  due  to  other  causes.  The  information  gained  can  be  used 
in  planning  an  experiment;  i.e.,  in  determining  specific  activity  a  radioisotope 
must  yield  a  desired  detection  limit  for  a  specified  counting  time;  the  infor¬ 
mation  may  also  be  used  to  determine  the  statistical  reliability  of  experimental 
data.  Established  methods  of  statistical  analysis  can  be  applied  in  the  estima¬ 
tion  of  the  error  due  to  random  variations  in  a  disintegration  rate  determina¬ 
tion  (R-l|  2,  3).  Only  a  brief  summary  of  the  most  basic  considerations  is  pre¬ 
sented  here. 

Use  of  Standard  Deviation 

It  can  be  shown  that  the  probability  of  occurrence  of  various  counts  which 

differ  from  the  true  average  rate  can  be  predicted  by  either  a  Poisson  or  a 

Gaussian  distribution  (by  the  use  of  approximations  which  are  normally  valid 

in  counting  radioactive  disintegrations).  The  standard  deviation,  o,  of  a 

single  measured  value  of  a  total  number  of  counts,  N,  may  be  expressed  as 
1/2 

a  «  (N)  .  This  is  valid  for  either  a  Poisson  or  Gaussian  distribution  when 

N  is  large.  The  standard  deviation  is  an  important  statistical  quantity  and 
possible  errors  corresponding  to  various  probabilities  may  be  estimated  by  fac¬ 
tors  of  it.  Thus,  the  uncertainty  of  a  counting  determination  which  is  associ¬ 
ated  with  a  particular  probability  of  occurrence  can  be  expressed  as  ±  K 
As  K  assumes  larger  values,  the  probability  that  an  error  will  exceed  ±  K 
diminishes.  The  values  of  K  which  correspond  to  a  number  of  commonly  used 
errors  are  shown  in  the  following  table. 
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Table 

Statistical  Errors  in  Counting  Determinations 

Probability  of  Observing 


Name  of  Error 

K 

Error  within 

Probable  error 

0.674 

0.500 

Standard  deviation 

1.000 

0.683 

Nine-tenths  error 

1.645 

0.900 

Ninety-five  hundredths  error 

1.960 

0.950 

For  example*  if  100  counts  are  recorded  in  a  given  time*  the  results  may 
be  expressed  as  N  t  0.674"!^  or  100  ±  6.7  counts.  There  are  fifty  chances  out 
of  a  hundred  that  the  error  in  the  number  of  counts  would  fall  within  a  ±  6.7 
counts  since  a  K  «  0.674  corresponds  to  a  0.50  probability  level.  It  is  empha¬ 
sized  that  the  uncertainty  refers  only  to  that  which  might  be  present  because 
of  random  fluctuations  in  the  decay  rate. 


Estimation  of  Error  in  Counting  Rate 

Usually  the  result  of  a  counting  determination  is  expressed  as  a  certain 
number  of  net  counts  per  minute.  This  value  can  in  turn  be  related  to  the 
number  of  disintegrations  per  minute  and  ultimately  to  a  concentration  value 
for  the  element  in  question.  The  net  counting  rate  is  simply  the  difference 
between  the  total  count  rate  and  the  background  count  rate. 

The  background  has  its  origin  in  various  natural  sources  of  radiation 
including  cosmic-radiation  and  the  naturally-occurring  radioactive  elements 
present  in  the  experimental  apparatus.  Since  a  low  background  rate  is  desirable* 
shielding  and  other  methods  are  used  to  minimize  it. 

The  standard  deviation  of  a  net  counting  rate*  which  is  a  difference*  is 
equal  to  the  square  root  of  the  sum  of  the  squares  of  the  standard  deviations 
of  the  total  rate  and  the  background  rate  (R-3) : 


(1) 
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Sine*  the  expression  for  the  standard  deviation  of  the  number  of  counts  is 
a,  ^ 

and  (where  t  is  counting  time), 

K  t 


I  R 

then  (where  R  is  counting  rate). 

Thus,  Equation  (1)  may  be  written  as  follows 


•  (-5-  •  -i-) 


1/2 


(2) 


Typical  Problem 

A  common  problem  which  arises  in  the  measurement  of  a  disintegration  rate 
is  the  following; 

The  count  of  a  sample  (including  the  background)  for  10  minutes  yielded  a 
total  of  426  counts.  The  background  was  counted  alone  for  10  minutes  and 
yielded  334  counts.  What  uncertainty  should  be  assigned  to  the  net  counting 
rate  of  the  sample  due  to  random  variations? 


Solution; 

The  net  counting  rate  of  the  sample  is 
426  ^^4 

Rg  ■  -yy  -  ■  9.2  counts/minute. 


The  standard  deviation  of  the  net  counting  rate  is  given  by  Equation  2. 


a 


R 


S 


/  42.6 

V  10 


10  / 


2.8  counts/minute. 


Thus  the  results  of  the  determination  may  be  expressed  as  9.2  t  2.8  counts 
per  minute  with  a  confidence  level  of  68.3  per  cent.  If  a  higher  confidence 
level  is  desired,  the  standard  deviation  should  be  multiplied  by  the  appropriate 
K  as  shown  in  the  Table. 


Jarrett  (R-2)  presents  an  analysis  designed  to  minimize  the  error  Intro¬ 
duced  by  the  background  by  use  of  an  efficient  distribution  between  the  total 
counting  time  and  the  background  counting  time.  Also,  Jarrett  presents  a  num¬ 
ber  of  curves  and  nomographs  which  are  convenient  to  use  in  estimating  the 
statistical  error  inherent  in  the  measurement  of  radioactivity. 
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DIFFUSION  THROUGH  SIO.  FILM  -  ANALYTIC  SOLUTIOI 


APPENDIX  B 


Mg  Diffusion  through  Si02  Film  -  Analytic  Solution 


The  diffusion  of  magnesium  in  gallium  arsenide  through  the  silicon  dioxide 
film  presents  the  problem  of  finding  the  concentration  at  the  interface  of  Si02 
and  GaAs.  This  problem  has  been  solved  analytically. 

Consider  the  SIO2  -  GeAs  semi-infinite  system  as  shown  in  Fig.  1. 


Fig.  1 

We  know  magnesium  diffusion  follows  the  complementary  error  function  distribution 
in  gall  ium  arsenide  and  therefore  concentration  -  independent  diffusion  coefficient 
in  GaAs  has  been  assumed, 

ThuSg  suppose  in  the  semi-infinite  region  -a  <  x  <  «  »  the  diffusion  coef¬ 
ficient  is  D|  in  the  region  -a  <  x  <  o  and  D2  in  the  region  o  <  x  <  »  ,  During 
diffusion!  two  operations  take  place: 

(1)  Inward  diffusion  in  the  region  -a  <  x  <  «  starting  from  the  surface 
x  ■  -a»  the  surface  of  silicon  diodeg  and 

(2)  Out-diffusion  from  the  region  o  <  x  <  ^  through  the  interface  x  >  o. 

The  solution  of  these  two  cases  should  then  be  added  electrically  to  give  the 
concentration  at  the  interface  Si02  *  GaAs, 
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The  boundary  conditions  for  these  two  cares  follow: 
Case  1 


(x,t)  -  (x,t) 

( -a  <  X  <  o,  t  >  o) 

(la) 

(x,t)  .  (x,t) 

(x  >  0,  t  >  o) 

(lb) 

o 

X 

o 

1 

o 

(-a  <  X  <  o) 

(Ic) 

C  {x,o)  ■  0 

(x  >  o) 

(Id) 

C  (-a,t)  - 

(t  >  o) 

(le) 

Llm  C  (x,t)  -  0 

(If) 

X  —a- 00 


D^Cjj(-o,t)  -  DgCj^C+o.t)  (Ig) 

C  (-o,t)  -  (1-H)  C  (♦o,t)  +  (Hi) 


Note:  Boundary  condition  (1h)  takes  into  account  the  difference  of  solid  solu 
bilities  of  the  diffusant  in  Si02  and  GaAst  the  concentration  in  GaAs  may 
approach  a  maximum  solid  solubility,  which  it  cannot  exceed.  H  is  the  step 

ft 

function. 

Case  2 


(x>t)  -  (x,t) 

(-a  <  x<  0{,t>  o) 

(2a) 

(*#t)  -  DgCjyj  (x,t) 

(x  >  0,  t  >  o) 

(2b) 

C  (x,o)  ■  0 

( -a  <  X  <  o) 

(2c) 

C  (x,o)  - 

(x  >  o) 

(2d) 

R.  C.  Wackwitz,  "Analytic  Solution  of  the  Multiple  Diffusion  Problems" 
Journal  Applied  Physics,  Vol.  33,  #9,  Sept.  1962. 
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vlicre  It  tbei  eoneentratlon  of  donors  in  GaAs. 


Lla  C  (x,t)  *•  (t  >  o) 

X  — *>oo 

(2e) 

(-*,  t)  -  o 

(2f) 

®1 

(2g) 

Dg  (+o,t)  -  K  (+o,t)  -C  (-o,t)J 

(2h) 

Hotoi  Boundary  condition  (2h)  assvunas  that  tbs  out-diffusion  from  QaAs 
to  SiOg  It  proportional  to  the  difference  between  GeiAs  and  SiOg  instan- 
taneout  surface  concentrations. 

For  small  ‘values  of  t,  'the  solution  of  equation  (l)  is 


“o.*.  <*•*)  -  ^  '“=0.A.  *  “ou.  (^) 


and  the  solution  of  equation  (2)  is 
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o  2 

/  Ua4-x  \  d  d,  Ai  V  ov2  2\  2ah  -hx  h  y 

erfc(^*!|~  +  hyj  +  -g  eric  “g  (l-^-2h  y  )  e  e  e  ' 

/2a -x\  ^  _  /  4a+x\2 

*  ^  °  ^  T  ^1"'° 

. A‘  •'''  *  “>')}  • 

'(-ij)  -  •“ 

:^)  ■  I  •"«  L'^)] 

-  2h^y^j 


erfc  +  ^yj+  I  (^  ■*■  J  ^ 

j.-lfp)"  *  .'^‘^^1 .  (i.hw.3rt-aV 

y  fiT 

^  «rfc  ■*■  hy^ 


h^x  >h 

Final  solution^  therefore |  looks  as  follows; 
•  C^  erfc 


®810, 


(x,t)  -  C^ 


(ig)«c.a.rfo(J^)-0.a».. 


\  cjT  /  o  \  «-jr  /  ■  \  f  • 

2y)-Ut.  »  J'"'!  2y}'L  ■  < 

X  ^'^^^GaAbI  ^  I  Ua+x\  “^x  r d  ii  ™.2. 

-  *  — —  \-wr  IS-  \  7  = 


,  /a.x\  [“^“x  ^  /W*x\  ""““x  r 

••"'  brj-  [;;;?-  *  — v-J'"'  tt  [ 


a' 

‘V) 


-hx^h^y^  erfc^^ 
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/2a-x\^  <•  “(^^) 

+  1  V;Uy;  .1  (2a.x)e 

(h-A?'  h^yAvL  3  hVV"  ^ 

♦  (W.)  .'  .'“  .  * 


A,(».‘) 


.r.c  {^). 


iia-fx  V  ^  1 

:w 

HC, 


(» 


3  a  erfc 

B 

+  _ 

E 


(¥)-  “ 

erfc/ 


w.  '"'("I;)' "  “o***  (^) 

'0..A.  “  (^1  ■  'x  ♦  r  •'''(-ij) 

|)  .rfc(-^  *  hjr).  g  |.rfc  (K^)-  (l-  a)erfc  (l!^) 

■  -  2»V).rfe(K^  .ay)  -  ^  ^ 

(l.hw.3ia..21.V)  .3*“  (“i^  *  '’»)]  <3b) 

“  '/^  ’  ^  ■  V°iV  •  “  ■  ^ 


‘(1*1  -h„x. 


a .  knl  ^  H 

l+(iH 

l-^H  ,  d  -  I  (1-m),  li  -  f  (1+^) 

Ul 


A.n^ 


Knowing  and  D^*  this  solution  will  give  the  concentration  value  at  the  inter¬ 


face  Si02  -  GaAs. 


This  solution  is  very  general  in  that  it  can  be  applied  to  any  semi-infinite 
composite  system  which  has  a  skin  or  surface  layer  of  finite  thickness.  Both 
media  should  have  concentration-independent  diffusion  coefficient. 

Further^  the  solution  given  by  expressions  (3a)  and  3b)  can  be  proven  to 
be  the  unique  solution  to  this  boundary  value  problem. 
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Xt  is  Interesting  to  note  th|t  the  values  of  f  (y>z)  where 

f  (Xj*)  ■  y  +  e^  e*  erfc  (y+z) 

”  * 
for  different  values  of  y  and  z  are  available  In  literature  . 

The  above  solutions  are  given  for  small  values  of  t.  The  general 
solutions  are,  however,  desired  when  the  diffusion  Is  carried  for 
greater  lengths  of  time. 

For  all  values  of  t,  the  solutions  can  be  written  as  follows: 


®^GaA8  **  ^ _ _  (  2na  +  x  +  2bl\ 

>  — ST- Cry— ;  •  “cT- '"n - TjF - ) 

Zn(x+a)  -Z„(x+a) 

+  e 

V 

be  °  +  ge 


xc. 


_  z„(x+a)  -Z  (x+a 


and. 


b  e^a*^  +  g  e'^a* 

JJ 

K(u+1)  (Z^a+2)' 

Where,  b  -  ^  (Z^^a  ♦  3) 

*  V«  'V'«  - 

°l^n  . 

K(u-l)  (Zna-2)~ 

®  ■  2K 

U'  *  -5)  J 

*  Stalts  and  Killer,  "Rate  Limitation  ?t  the  Surface  for  Impurity  Diffusion 
in  Semieonduetors" ,  Bell  Telephone  System  Technical  Publications  Moiu)graph  3^^9> 
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ftndf  Zjj  ( B  ■  2^  3i  •  •  •  • )  ftr#  the  zoots  of  the  equation 
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